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ABSTRACT 


Solutions are derived for adhesive-bonded joints of non-classical geometries. 
Particular attention is given to bonded doublers and to selective reinforcement 
by unidirectional composites. Non-dimensional ized charts are presented for the 
efficiency limit imposed on the skin as the result of the eccentricity in the 
load path through the doubler. It is desirable to employ a relatively large 
doubler to minimize the effective eccentricity in the load path. The transfer 
stresses associated with selective reinforcement of metal structures by 
advanced composites are analyzed. Reinforcement of bolt holes in composites by 
bonded metal doublers is covered quantitatively. Also included is the adhesive 
joint analysis for shear flow in a multi-cell torque box, in which the bond on 
one angle becomes more critical sooner than those on the others, thereby 
restricting the strength to less than the total of each maximum strength when 
acting alone. Adhesive plasticity and adherend stiffness and thermal imbalanc- 
es are included. A simple analysis/design technique of solution in terms of 
upper and lower bounds on an all-plastic adhesive analysis is introduced. This 
is far simpler than the more precise elastic-plastic adhesive analysis and, in 
most cases is of adequate accuracy. An analysis of tapered-lap bonded joints 
is included and this shows how to alleviate the peel-stress problem character- 
istic of thick uniform double-lap joints. The tapered-lap joints remain fully 
efficient for thicknesses above which the uniform double-lap joint becomes 
inefficient. Illustrative examples are included throughout the text. 
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SUMMARY 


It is not usual to find adhesive bonded joints in aerospace structures which 
conform precisely to the geometry and load conditions of the classical classes 
of joints on which tests are performed. The objective of this report is to 
apply to practice some of the analysis/design capability developed by the 
elastic-plastic formulation of the classical (double-lap, single-lap, stepped- 
lap, and scarf) joints. Particular attention is devoted to bonded doublers, 
with a range of end support conditions. Non-dimensional ized charts are present- 
ed for the efficiency limit imposed on the skin as the result of the eccentric- 
ity in the load path through the doubler. It is desirable to employ a relative- 
ly large doubler to minimize the effective eccentricity in the load path. In 
addition, consideration is given to the selective reinforcement of metal struct- 
ures by unidirectional composites and of composites by bonded metal doublers 
around bolt holes. In the former case, the key issue is the load transfer at 
the end of the composite reinforcement while, in the latter case, the usual 
problem is that of load transfer between the metal doubler and the composite 
part. Also to be found in this report is an analysis of the bond stress and 
strain distributions for shear flow in a multi-cell torque. A technique of 
alleviating the peel-stress problem for thick uniform double-lap joints is des- 
cribed and analyzed. The optimum tapering of the outer adherends is shown to 
eliminate the peel problem while adding 24 per cent to the bond shear strength. 

Adhesive plasticity and adherend stiffness and thermal imbalances are accounted 
for. A simple analysis/design technique of solution in terms of upper and lower 
bounds on an all-plastic adhesive analysis is introduced. This is far simpler 
than the more precise elastic-plastic adhesive analysis and, in most cases, is 
of adequate accuracy. Indeed, in many cases such solutions are simpler than a 
perfectly-el asti c method. Consequently, analytical solutions can be obtained 
for a far graeter range of joint configurations and complex load conditions than 
is possible with more precise methods. 

Illustrative examples are included throughout the text to explain how to employ 
the analyses derived. 
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1. INTRODUCTION 


It is unusual in the design of aircraft and space vehicles to encounter an 
adhesive-bonded joint of precisely one of the classical families of analysis 
such as the double-lap joint, the single-lap joint, the scarf joint, and the 
stepped-lap joint. It is even more unusual to encounter simultaneously a pure 
form of load application, such as tensile lap shear which is the basic test 
condition. The differences in load capacities for compressive lap-shear and 
in-plane shear are discussed in companion reports (References 1,2, and 3). 

The purpose of this report is to apply the basic analysis techniques established 
in References 1, 2, and 3 for double-lap single-lap, and scarf and stepped-lap 
joints, respectively, to joint configurations more frequently encountered in 
aerospace practice. 

Perhaps the most numerous applications of adhesive bonding are to be found in 
the form of bonded edge doublers for flush mechanical attachments. Such applic- 
ations share many of the governing equations for unsupported single-lap joints. 
The limiting characteristic is usually not that of adhesive shear or peel but of 
the non-uniform load distribution across the thickness of the main sheet just 
outside of the doubler. It has long been known that the doubler must be flex- 
ible in bending at its tip to diffuse load transfer gradually. This is why 
feathered edges and fingered doublers are so widely employed. The analyses here 
confirm the need for such procedures and provide a rational basis for design. 

Another application, of growing importance in the form of metal structures 
selectively reinforced with advanced filamentary composites, is that of a 
bonded-on area doubler, usually unidirectional unless it is thick. Such 
doublers are quite distinct from edge doublers inasmuch as the latter serve 
principally to reduce a stress level locally while the former are more of the 
principal load path type of structure. As such it is usual to minimize the 
eccentricities in the load path, to provide moment-resistant supports, and aim 
for the much higher efficiencies of double-lap joints. The key problem with 
this class of bonded structure is usually that of load transfer at the ends. 

If the area of the bonded-on doubler is too large in proportion to the bond 
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area, the weak link will be the shear capacity of the adhesive at the ends of 
the doubler. 

A common problem in composite structures is that of bonded metal doublers in an 
area of mechanical attachments. An analysis technique and illustrative example 
are provided. 

The fourth class of non-classical adhesive-bonded joints is concerned with the 
transfer of in-plane shear loads. The pure form of this problem, from one 
member to another, is covered in References 1 and 3. (The single-lap case is 
governed by the analyses of Reference 1 for this situation.) The analytically 
more challenging problem arises in multi -cell torque boxes, such as aircraft 
wings. The best known illustrative example is that of a wing skin to spar web 
attachment, with a change in skin shear flow at the intersection. The boundary 
conditions for each leg of such a joint differ from those of a single-transfer 
joint, but the governing equations are of the same form. A related problem also 
covered by the analysis in this report (with appropriate boundary conditions) 
is that of the shear flow between a spar (beam) web and a concentrated cap. 

Approximate analysis techniques for combined loading on joints are suggested. 
These are based on the concept of the maximum adhesive shear strain being the 
square rooot of the sum of the squares of orthogonal strain components. 

In the more complicated joint configurations it does not prove practical to 
perform precise elastic-plastic analyses. It is quite effective to obtain 
simpler solutions in terms of upper and lower bounds by means of perfectly- 
plastic adhesive analyses. The simplification arises from the elimination of 
the "boundary" conditions at each of the elastic-to-plastic transitions. This 
simplification is even more powerful than just a reduction in number of boundary 
conditions because the location of such transitions is not known at the start of 
an analysis and, furthermore, it shifts as the load intensity is changed. (See 
References 1, 2, and 3). Further justification for this technique is provided 
by the fact that, for ultimate static load, the upper and lower bounds are 
quite close together for the ductile adhesives used in subsonic commercial 
aircraft structures. Even though such bounds are further apart for the brittle 
adhesives needed for high- temperature zones on supersonic aircraft, they 
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represent a considerable improvement over the blanket application of a uniform 
shear stress developed from single-lap shear tests. This report, in conjunction 
with the companion References 1, 2, and 3 provides all the background and 
techniques necessary for a more precise elastic-plastic analysis if ever needed. 

The report concludes with a description of a powerful technique to obtain peel- 
stress relief for thick uniform double-lap adherends. This consists of tapering 
the outer adherends to make them flexible at the tip in order to prevent the 
peel stresses from ever developing. This technique has been used in practice 
for a considerable time. By means of the fully-plastic analysis method 
outlined above, it was possible to determine the optimum proportions to not 
only achieve peel -stress relief but also to simultaneously maximize the shear 
transfer capacity of the tapered-lap bonded joint. 
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2. ADHESIVE-BONDED DOUBLERS 


One of the most frequent applications of adhesive bonding in aircraft struct- 
ures is that of bonded edge doublers. Two important applications are the build 
up in thickness to permit the use of flush fasteners in thin gage sheet metal 
and the fingered or feathered-edge doublers around panels subjected to accoustic 
fatigue loads. In either case, the bond itself is not usually the limiting 
element in the strength of the structure. The adhesive shear load can usually 
be developed in about half an inch or so of plastic zones (while the overlaps 
are typically very much greater) and typical design practice excludes substant- 
ial differences in thickness at the edge of the doubler (and thereby minimizes 
any potential peel-stress problems). The dominant limiting feature is the 
eccentricity in load path when a doubler can be fitted on one side only of the 
structure as is the case for external skins. The relevant analyses are akin to 
those for the single-lap joint (Reference 2). 
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2.1 Adherend Stresses in Adhesive-Bonded Doublers 


It is customary to make fingered doublers less thick than the skin to which 
they are bonded (typically 0.6 times as thick or one skin gage less if very 
thin). Therefore, with reference to the single-lap joint analysis of Reference 
2 , equation (152) is of interest here rather than the moment at the other end 
of the single-lap joint which dominates that analysis. The distance (c - e) in 
equation (152) of Reference 2 refers to the distance between the point of 
inflexion (m = o) and the edge of the doubler and may consequently be less than 
the total overlap. Here, the end-support condition selected for analysis is 
that of the skin/doubler combination being built in. Also, the analysis of 
Reference 2 is adapted to provide for simply-supported ends. A comparison of 
the two sets of answers so derived indicates that the former is related to the 
latter by a factor of 2 (to some five or six significant figures) in the extent 
of the overlap. Other end-support conditions can therefore be analyzed by 
interpolation. 


The geometry governing the present analysis is defined in Figure 1. The load 
is assumed to be reacted at the neutral axis of the skin/doubler combination. 

It is easily shown that this stiffness "centroid" is located at a distance from 
the middle of adherend 4 of 


EC 4 


|(ti + t4 + 2n) 



( 1 ) 


Throughout the adherend 4, outside the overlap, the longitudinal stress result- 
ant (force per unit lateral width) is uniform, at the value P of the applied 
load. The stress couple M 4 , per unit width, is then defined by the equilibrium 
equation 


M u = P 


(ti + + 2n) /X4 


( £4 + 1 ) ( 1 + 


E4t 4 V 

E^t 


r 


0 - 


W4 


for -£4 i X4 £ 0. 


( 2 ) 


The classical theory for the infinitesimal deformation of thin, cylindrically 
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bent plates then yields 


d 2 W4 M4 

p 

(tj + t 4 + 2n) /x 4 \ 

dx4 2 D4 

d 4 

/ E 4 t 4 \\ 2 / 

(£4 + £)|l + 

\ Ejt ! / 


whence 


w 4 = A4 cosh( £4X4 ) + b 4 sinh( ^4X4 ) + 


(ti + t 4 + 2n) /x 4 




(1 ♦ ^ 

\ E x t x / 


s- 


( 3 ) 


(4) 


where 


P 12P( 1 - v 4 2 ) , % 

u 2 = r ■ < 5 > 

D4 ( ) 4Ei + 't 4 

The assumption of negligible moment restraint at X 4 = 0 (or, equivalently, that 
£4 is so large in comparison with £ d that the precise nature of the end condi- 
tion is immaterial) prescribes that 


a 4 = 0 


( 6 ) 


The behavior of regions 2 and 3 (Figure 1) is governed by precisely the same 
equations as for the single-lap joint analysis (Reference 2), with the excep- 
tion that, here, at s = 0 the boundary conditions differ. The equations are 
set up simply as 


dM 2 


ds 

dM 3 


ds 



0 


0 



dT 2 

+ t = 0 

ds 

dT 3 

- t = 0 

ds 


(7) 


(8) 
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dV 2 

°c = ° ^ 


ds 



dV 3 

a=0 ) 


ds 

U 

d 2 W2 

m 2 

12M 2 (l — v 3 2 ) 

ds 2 

Dl 

(k b )iEi t i 3 

d 2 w 3 

m 3 

12M 3 (1 - v 4 2 ) 

ds 2 

d 4 

(k^) 4 E 4 t 4 3 


As in Reference 2 it is necessary to use the approximations 

dv2 dW3 

V2 " W 3 , = 

ds ds 


( 9 ) 
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(ID 


while maintaining the distinction between the higher derivatives. From the 
equations above. 
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(12) 

(13) 

(14) 


(15) 


and, since for an adhesive obeying a linear elastic law in tension (peeling) 


d 4 (w 2 - w 3 ) 
ds* 4 


a c ( w 3 - W2) 



IDj D 4 ' n 


1 \ E c / 

— nkw2 - w 3 ) + ( 
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2Dj 
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(16) 

(17) 


Equations (14) and (17) indicate coupling between the shear and peel stresses 
in the adhesive. In order to obtain an explicit solution expeditiously, an 
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approximate solution will be adopted with sufficient integration constants to 
satisfy the dominant boundary conditions. The simplest such solution, used 
also in Reference 2, is 

^-(w 2 + W 3 ) - A 23 S 3 + B 23 S 2 + C 23 S + F 23 . (18) 

The conditions at s = 0 for a built-in end of the skin/doubler combination 
permit the setting of 

C 2 3 = F 23 = 0 . (19) 

(For simply-supported ends, the corresponding relations would be 

B 2 3 F 23 = 0 . ) ( 20 ) 

As in Reference 2, it is assumed that 


dw 2 dw 3 


ds ds 


at s = 0 and s = a. 


but that the distinction between the second derivatives is of paramount import- 


ance since 


d 2 W3 M 


ds 2 Dl 


= 0 at s = n. 


(Note the sign convention for positive values of m q , the critical moment in the 
adherend.) The pertinent boundary conditions, at s = & d> X 4 = -£4 are 


W4 = - B4sinh( £4^4 ) - 


ti + t4 + 2n 


£4 + l 


d J 1 + 


= i(w 2 + w 3 ) = A 23 (0 3 + B 23 U,) 2 , 


= B 4 ^4COsh(5L f £4 ) 


+ t 4 + 2 n 1 1 


^4 + J(i 


E 4 t 4 \ 2 


d(w2 + W3) 


= 3A 23 (^ d ) 2 + 2B 23 U d ) , 
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d 2 w 4 M 

= -2. = _ B 4 C 4 2 sinhU 4 4 4 ) 

dx 4 2 D 4 


/d 2 (w 2 + w 3 ) d 2 (w 2 - w 3 ) 


)■ 


M 

2D 4 


\ ds z ds^ 

and, for sufficiently large £; 4 £ 4 , 

sinh( 54 ^ 4 ) - cosh ( £ 4 ^ 4 ) - ie^ 4 S ^ 4 ^ 

2 

By a process of elimination it is established that 
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/ t 2n\ 

1 + - + ~ 
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(25) 


(26) 


(27) 


(28) 


Equation (27) holds for built-in ends. It follows from equation (152) of 
Reference 2 that, for simply-supported ends. 


M 

o 


1 + (c A ) + i(S 4.) : 
s d 6 s d 



2n“ 


1 + 


E t 
s s 

Vd 


(29) 


2 

In these equations, the subscripts d and s refer, respectively, to the doubler 
and skin, being equivalent to the subscripts l and k. The bending stiffness 
parameter £ s for the skin is given above as 

P 12P(l - v 2 ) 

? 2 =— = — (31) 

S D (k. ) E t 3 

S D S S S 

In this, the coefficient k serves to uncouple the bending and extensional 
stiffnesses so that the analysis covers filamentary composite materials as well 
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as the isotropic metals. (In the computer programs listed in the Appendix, the 

ratio E /E is set equal to unity. That is, the skin and doubler materials are 

s d 

the same. This restriction is in accord with common practice, but can easily 
be over-ridden in a specific case. Likewise, the assumption in equation (28) 
that the doubler is small in comparison with the extent of the skin can also 
be refined if desired.) The maximum stress in the skin, at the edge of the 
doubler, follows from equation (27) or (29) as 

P 6m 

a =— + -^=a (1 + 3k) . (32) 

max ^ ^2 av S 

s s 

Equations (28), (30) and (32) are programmed to be solved by iteration in the 
digital computer program listed in the Appendix. Representative solutions, in 
non-dimensional ized form, are shown in Figures 2 and 3. It is apparent that, 
for all practical purposes, the point of inflexion in the built-in case is in 
the middle of the doubler since doubling the i/t ratio with respect to the 
simply-supported case yields the same joint stresses when computed for the 
built-in case. It can be seen clearly that the effective length of the overlap 
should be as large as practical. Trying to save weight by skimping on the 
doubler has the effect of reducing significantly the allowable skin stress 
(a /a < l) so that a major weight penalty can be suffered for the entire 

avg max _ . 

panel. This is why fingered doublers are customarily applied with quite long 
overlaps rather than shallow scallops. The analysis above governs the fatigue 
failure mode illustrated in Figure 4. The cracks originate on the doubler side 
of the skin at the extremities of the fingers. Such cracks do not tend to 
occur on the (sloping) sides of the fingers because no clearly definable axis 
for bending of the skin exists there. The experimental evidence is strong that 
doublers should be either fingered or have feathered edges to prevent the form- 
ation of a long continuous crack in the skin adjacent to the edge of the 
doubler. 

It will be noticed that nowhere in the analysis above is there any involvement 
of the adhesive properties. Consequently the solution can be applied also to 
integral doublers formed by chem-milling of a thick skin. In doing so, 
however, it must be recognized that there is an additional stress concentration 
due to the abrupt change in thickness which is not eliminated for the integral 
doubler as it is for the bonded doubler by the layer of adhesive. 
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EXAMPLE 1: 


To illustrate the use of the analysis above, consider a 2024-T3 aluminum skin 
0.025 inch thick with a doubler 0.015 inch thick. The required skin efficiency 
is 90 per cent, so the minimum overlap is read from Figure 2 as 

l = 0.025 x lL.85 / 0.22168 = 1.67 inch. 

(An increase to z = 2.65 inch overlap would be needed to raise the skin 
efficiency to 95 per cent.) 
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2.2 Adhesive Shear Stresses in Adhesive-Bonded Doublers 


Consideration of adhesive stresses is not necessary in most bonded (edge) 
doubler applications because the analysis above for skin efficiency tends to 
over-ride such considerations. Nevertheless, the bonded doubler serves to 
illustrate simply the use of modified all-plastic design concepts for 
structural bonded joints. 

The natural desire for high adherend efficiency is satisfied by minimizing the 
effect of any eccentricities in load paths. The accomplishment of this mini- 
mizes the differences between the actual adhesive stress distributions and 
those of a supported (double-lap) joint. Therefore, instead of adapting the 
more complex analysis of Reference 2 to the boundary conditions of bonded 
doublers, the simpler analysis techniques of Reference 1 are employed. The 
analysis presented here (see Reference 4) accounts for adhesive plasticity, 
while assuming that the skin and doubler remain elastic. An algebraic solution 
is derived which is sufficiently simple for design use. The solution describes 
the adhesive shear stress distribution between the skin and doubler, identify- 
ing the critical region (a narrow strip along the edge of the doubler) and 
those other regions in which small manufacturing defects can be tolerated with- 
out impairing the structural efficiency of the design. This analysis shows how 
to size minimum bond areas by estimating the extent of plastic adhesive zone 
necessary to transfer the load (and check that the adhesive is capable of gen- 
erating such strengths for the given adherends) and adding to this the distance 
required to build up the elastic adhesive stresses from zero to the plastic 
level. This approach is more meaningful than that of sizing in terms of a uni- 
form bond shear stress, no matter how the "allowable" was arrived at. 

Peel stresses in the adhesive (or interlaminar tension stresses in an adjacent 
composite adherend) become a progressively more severe problem for either short 
overlaps or thick sections. Usually these conditions are not encountered in 
practical doubler design. In the event that they are, the single-lap joint 
analysis of Reference 2 can be modified easily in much the same manner as for 
the adherend analysis in Section 2.1 above. 
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Figure 5 depicts the geometry and nomenclature for the analysis of an arbitrary 
adhesive-bonded skin/doubler combination. A uniform tensile load P (per unit 
width) is applied to the skin and taken out at a row of close-spaced rivets 
adjacent to the outer edge of the skin/doubler combination. A compressive load 
would be equally well accounted for by this analysis, but it is necessary here 
to establish a sign convention. In order to define the problem sufficiently to 
permit a solution to be obtained it is assumed here that, at the row of rivets, 
both skin and doubler act is unison, each transferring load to the rivets in 
proportion to their extensional stiffness. It is also assumed that the doubler 
extends sufficiently beyond the row of rivets for the local stress variations 
around the rivets not to affect the analysis in the region of adhesive shear 
load transfer. The row of rivets are replaced mathematically by a thick 
extension beyond the end of the original skin and the doubler, as depicted in 
Figure 5. 


The conditions for horizontal force-equilibrium for a differential element dx 
within the adhesive bonded area are 


dT dT 

+ x = 0 , — — - x = 0 , (33) 

dx dx 

in which the subscript s refers to the skin and d to the doubler. The stress- 
strain relations for the skin and doubler yield 


d6 


T 


d6 . 


¥d 


dx E t dx 

s s 

The adhesive shear strain y is approximated by the relation 


Y = 


(6 - 6 ) 
s d 


(34) 


(35) 


and, within the elastic region (of unknown length s), the adhesive shear stress 
is assumed to be 


G(6 - 6 ) 

T = Gy = — - — — = f(x) (36) 

n 

while, throughout the remaining plastic region (of unknown length w), the 
adhesive shear stress is assumed to be constant, at the plastic adhesive shear 
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stress 


t = r = constant 
p 


(37) 


To solve the governing differential equations, one eliminates <5 and 6 between 

S 0. 

equations (34) and (35) and then uses equations (33) to eliminate t and t , 
yielding the differential equations 


dY 

dx 


E t 
- s s 


VdJ 


(38) 


d 2 y 


E t 
■ s s 


E t 
d d 


T = 0 


(39) 


dx z n | 

Within the elastic region, which is assumed to extend to the rivet line because 
of the small relative displacement in that region, equation (39) becomes 


d 2 y 


dx z 


X 2 y = 0 


A 2 = 


n\E t 
s s 


E dv 


In the elastic solution, 

Y = A sinh(Ax) + B cosh(Ax) 


(40) 


(41) 


the constant B can be set identically equal to zero because of the assumed 
absence of relative displacement across the bond line at the row of rivets. 
Equation (38) then predicts that the loads in the skin and doubler are not 
proportional to their respective stiffnesses at the row of rivets. In other 
words, the shear deformation in the adhesive prevents the doubler from fully 
developing a load in proportion to the relative stiffness of the skin and doub- 
ler. Provided that the doubler is not excessively narrow, this effect is 
inconsequential, as is confirmed by the analysis below. 


In the plastic region, of assumed length w at the inner end of the doubler, the 
solution of equation (39) is 

Y = T p S 3 + C? + F (42) 

in which dx = d£ , the origin for 5 being at x = +s. 


The constants B, c and f and the unknown s or w are found by satisfying the 
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boundary conditions 


Y = 

Y e 1 

at 

X = s, £ = 0, 

(43) 

Y = 

Y e + 

Y P 

at £ = v, 

(44) 

. 

dx 

- ix 

dC 

at 

x = s, ? = 0, 

(45) 

. 

d? 

P 

“ E t 

c < 

T) 

c 

3- 

II 

(46) 


S S 


of which equation (45) ensures continuity in the adherend stresses. Hence 


F = Y = T /G , 
e p * 

A = y /sinh(As) = t /[G sinh(As)] , 
6 P 


Xt 


C = 


G tanh(As) 

From equation (45) it now follows that 


Aw + 


tanh( As ) 


1 ^ 

+ 2 - 2 - 


tanh 2 (As) Y 


so that, from equation (46) 


V 1 


P = — ^Aw + 


1 + 

tanh(As)/\ E 


t \ 

J ’ 


whence 


tanh 2 (Xs ) 


+ 2 — = 


XP 


/ E x \ 

T (l + 

_ p ' Vd/J 


(47) 

(48) 

(49) 


(50) 


(51) 


(52) 


Precise evaluation of these equations requires a digital computer program 
but, for sufficiently long overlaps, tanh(As) + l, whence 


[\v + l) 2 -> 1 + 2— •+ 


AP 


/ E t ' 
/-i s s 

T p( E,t. 


d d; 


(53) 
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The failure criterion for the adhesive is taken to be the exceeding of its 
total shear strain capacity (y e + Y p )> which has been shown to be a realistic 
approach in Reference 1 . 

Just as in the analysis of bonded lap joints (References 1 and 2), the e/tent 
of the plastic adhesive zone is limited by the adhesive and adherend properties 
and does not increase indefinitely with increasing applied load. ,.rst of the 
load transferred by the adhesive bond between the skin and the doubler passes 
through a single narrow band at the edge of the doubler. This same narrow band 
follows the contour of finger doublers, as shown in Figure 6. Manufacturing 
defects in the bond away from the edges do not impair the load transfer 
capability of the bond. 


Attention is now focused on the not quite ideal sharing of the load between the 
skin and doubler, as indicated above. The integral of the adhesive shear 
stress over the extent of both the elastic and plastic regions is 


S V 

/ T dx / 


= 


A sinh(As) 


-[cosh (As) - l] + t 


w 


(54) 


o o 

If the skin and doubler were to share the load, at the row of rivets, in 
proportion to their respective stiffnesses, the load per unit width in the 
doubler (transferred through the adhesive bond) would have to be 



which, from equation (51), is precisely equal to 


(55) 


A tanh(As) 


+ x w 
P 


(56) 


A comparison with equation (54) indicates that it is realistic to use the 
approximation (55) whenever cosh(As) »l which is practically always, because 
of the characteristically high values of A. 


The analysis above is now capable of re-interpretation in a form providing a 
meaningful comparison with the analysis of lap joints. Corresponding to equa- 
tion (50), the equivalent analysis of the general double-lap joint yields equa- 
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tion (A. 55) of Reference 1: 


X"b + tanh 



tanh 2 



(57) 


In lap joints it is the extent of this plastic adhesive zone which limits the 
ultimate strength that can be attained by the joint. Likewise, for bonded 
doublers, it limits the load that can be transferred from the skin to the 
doubler without failing the bond. In the latter situation, however, the poten- 
tial bond strength usually exceeds by far the strength of the adherends. At 
relatively low load intensities, the full plastic capability of the adhesive is 
not fully utilized. As the load is increased, a maximum load capacity is 
attained which fully extends the adhesive in shear. At this load level, no 
increase in the size of the doubler can increase the load transferred. Equa- 
tions (54) to (56) and (50) indicate that the load transferred through any bond 
of realistic length is equal to 


1 + 


E t 
s s 

E-t- 
d d 



cosh(As) - 1 


sinh( As ) 




(58) 


which is proportional to the square root of the strain energy of the adhesive 
in shear, just as for bonded lap joints. This explains why ductile adhesives 
are so much stronger than brittle ones despite the greater maximum shear stress 
of the latter. 


For practical stressing of the adhesive in shear, equation (58) contains all 
the pertinent information. The design sequence is as follows: 

( E t 
1 + S S 

Va 

(2) compare (1) with the maximum bond strength 



which is limited by the adhesive shear strain capacity (Y e + Y p ) and 
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( 3 ) if the comparison in ( 2 ) indicates that the adhesive is suitable, 

determine the minimum width of the doubler, from the row of rivets to 
the nearest edge of the fingers on the doubler as at least the sum of 
the widths of the plastic and elastic adhesive zones. For this pur- 
pose it would appear to be good design practice to base the load on 
the ultimate strength of the skin with an adequate margin of safety. 

A simple and generally not unduly conservative approach is to assume 
that tanh(Xs) - l for the elastic region, so that s = 3 /x, and that 
all the load is carried in the plastic region, so that 


W 


= P / 


1 + 


E t ^ 
s s 


E.t J 
d d' 


(4) Then, £ . . = s + v 

minimum 


EXAMPLE 2: 


Consider a 0.030 inch 7075-T6 aluminum alloy skin with a 0.020 inch doubler of 
the same material, bonded together with Epon 951 adhesive, for which t = 6 ooo 

ir 

psi , y p /y = 20 , g = 6.0 x io 4 psi , and n = 0.005 inch. For these materials. 


x = 



6.0 x IQ 4 

5 x 10- 3 


1 /10 3 
10 7 ' 30 


+ 


ion 

20 / 


10.0 . 


Taking a = 75 ksi for the aluminum (for which also E = 10 7 psi), the maximum 
bond load per inch is 75 x io 3 x 30 x 10 -3 / (l.O + 1 . 5 ) = 900 Ib/in. The 
adhesive bond capability is (6000 / 10 ) x Ai = 3842 lb/in., so the adhesive 
has more than adequate strength. 


The minimum width of the doubler (see Section A-A in Figure 6 ) then follows as 

3 . 900 . , _ inch 


£ = 


= 0 . 1+5 


10 .~0 6000 

The fingers can then be proportioned according to current accepted design 
practice. The minimum width £ above is likely to be so small in general that 
manufacturing tolerances and the adherend efficiency considerations above may 
well dictate a much greater length. 
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3. SELECTIVE REINFORCEMENT BY BONDED-ON DOUBLERS 


In the application of advanced composites to aerospace structures in the form 
of selective reinforcement, the objective is to minimize the total cost of the 
composite by restricting its application to those situations in which it is 
used most efficiently. These are the unidirectional reinforcement of beam 
caps, strut flanges, stringers and longerons. The key problem in such applic- 
ations is any load transfer from the composite to a metal end fitting. The 
brittleness of the composite has led to the acceptance of the concept of using 
only adhesive bonding between the composite and the end fitting and employing 
mechanical fasteners between the end fitting and the adjacent structure. The 
direct bolting of composite to the structure would necessarily require cross- 
plies to develop adequate bearing strength and, with the associated stress 
concentrations around the holes, this does not represent a particularly effic- 
ient application of advanced composites if the bolt holes run along the entire 
length of the part. A related concept, which has not received the attention 
it merits, is that of a bonded-on unidirectional reinforcement with cross-plies 
interleaved at the ends which can then be bolted as well as bonded. The effic- 
ient use of interleaved metal reinforcement or boron film reinforcement is 
limited to fairly thin sections because of non-uniform load sharing problems 
associated with the many rows of bolts needed to transfer load in to and out of 
a thick section. Under such constraints, it becomes apparent that reliance on 
pure adhesive bonding will continue in such applications. 

There are two mathematically distinct types of problems in this classification: 
the thin laminates with uniform thickness and the thicker ones with a scarfed 
(or stepped) end to effect the load transfer over a greater length. These are 
treated below in turn. Because of the interest in boron-infiltrated composite- 
reinforced structures (Reference 5), the solutions will be presented directly 
in terms of section areas per unit bond width rather than in terms of "equi- 
valent adherend thicknesses". 
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3.1 Uniform Thickness Doublers 


Starting with the uniform thickness adherends, the basic problem is that of how 
much composite can be bonded on before it splits off. This problem arises as 
the result of thermal stresses induced by the dissimilar coefficients of therm- 
al expansion for metals and composites and the fact that the operating temper- 
ature is almost invariably much less than the bonding (curing) temperature 
required for good environmental resistance. Once this limit is set, one can 
proceed to compute how much mechanical load can be transferred and how best to 
do so. The pertinent variables are depicted in Figure 7, in which the link 
supports are meant to indicate the lack of eccentricities and of any deflec- 
tions perpendicular to the x-axis. Precise elastic-plastic analysis will show 
that, to all intents and purposes, most of the adhesive load is developed in 
the end zones of the joint while the bulk of the structure is subjected to 
essentially uniform stresses. The distance d shown is independent of the total 
length for long overlaps but is characteristic of the cross section and the 
materials and the temperature differential at. The governing equations, for 
the sign convention employed in Figure 7, are: 

for equilibrium 

dTi 

- vt = 0 , (59) 

dx y 

and 

dT 2 

+ WT = 0 , (60) 

dx ^ 

for continuity 

y = ($i - 62 ) / n , (61 ) 

while, for the stress-strain relations for the adherend materials 


d6 : T X 

= + ai AT > (62) 

dx EjAi 

and 
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( 63 ) 


d .6 2 T2 

= + a 2 AT . 

dx E 2 A 2 

The governing boundary conditions for the adhesive are that 

t = t for 0 i x < d , 
p 


and 


t = 0 for x > d , 

Y = y at x = 0 , 

P 


(64) 

(65) 

( 66 ) 


Y = 0 at x = d . (67) 

The joint is symmetrical with respect to the centerline at x = L/2. 

In the range d £ x < L/2, there is no differential displacement across the 
adhesive layer so that, from equations (62) and (63), 

Ti T 2 

+ dj AT = + a 2 AT . (68) 


At the elastic-to-plastic transition at x = d, equilibrium requires that 

Tj = -t 2 . (69) 

The simultaneous solution of equations ( 68 ) and (69) leads to the results 


T 1 _ ( “2 “ a l ^ AT /(e 1 Ai + E 2 A 2 ) 


■“•2 


(70) 


This same result is deduced for the location x = l/2 by a perfectly elastic 
analysis, in association with an exponential increase in adhesive shear stresses 
towards the ends of the joint (see Reference 6 ). 


Turning now to the adhesive stresses, equations (59) and (60) require that the 
forces t x and t 2 be reduced linearly to zero at the free edges of the specimen 
by a uniform shear stress over the distance d. It follows that the length of 
plastic adhesive zone at each end of the overlap is given by 

d. | (.2 -.,)«! + M ; (71) 

WT^ / \EiAi E2A2/ 

in which the modulus signs are necessary to ensure that the length d is positive 
without prior knowledge of the relative magnitudes of 04 and a 2 . It is now 
appropriate to relate this distance d to the shear strain capacity of the 
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adhesive, which serves as the failure criterion. In the end zone, from equa- 
tions (59) through (63), 

Y = ifs, - i 2 ) - C + i (a, - d£ )ATx + • (72) 


At x = o. 


y = C = y e + Yp 

(x = 0) 

(73) 

while, at x = d. 



y = 0 

(x = d) 

(74) 

for the upper bound solution and 



Y = Y. 

(x = d) 

(75) 


for the lower bound solution. 


While the condition (75) is obviously associated with some displacement incom- 
patibility for x > a, equation (65) requires that there be no change from the 
forces in equation (70) for long and moderate overlaps. For the upper bound 
solution, then, equations (71) through (74) indicate that 

( Y e + Y p ) = ^ [ l { «2 - «1 )A T| ] , (76) 

whence the failure condition 

2 , p m( \ + y (ijAl + e^at) = [( “2 - a i ) at] 2 . (77) 

Various dominant influences are apparent in this relation. The thermal stress 
terms appear to the second power, so their influence becomes progressively more 
severe. With all other variables held constant, the left-hand side indicates 
a definite limit in adherend stiffnesses that can be loaded up, without failure 
of the adhesive, by a given adhesive and width of bond. The greater is w in 
proportion to E 1 A 1 and e 2 a 2 , the higher the temperature differential that can 
be withstood. The important role of adhesive ductility is again evident. The 
quantity nx_(Y + y ) represents the strain energy per unit bond area for the 
adhesive and this alone characterizes the influence of the adhesive. The same 
importance of the adhesive strain energy is evident in Reference 1 for double- 
lap joints. The lower-bound equivalent of equation (77) is 
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( 78 ) 


2t wr, Y (-V + fV) = I(«2 - «l)AT] 2 . 

P P\EiAi E 2 A 2 / 

Equations (77) and (78) are appropriate for the AVCO boron-infiltrated extrusion 
concept (Reference 5) because the circumferential symmetry prevents the devel- 
opment of peel stresses at the ends of the overlap. For the exposed type of 
reinforcement illustrated in Figure 7, it is probable that peel stresses or 
associated interlaminar tension stresses impose a more severe restriction than 
equations (77) or (78). This kind of problem is discussed in Section 6 of Ref- 
erence 1 and imposes a definite limit on the thickness of doubler which can be 
bonded on without having a scarfed end. For a uniform thickness of reinforce- 
ment, designated by the subscript o, this thickness is 


4 - 

°o 


E n 
o 


a cY 


3E '(1-v 2 )\t , 
c p 


(79) 


in which a is the maximum allowable peel or interlaminar tension stress, E ’ 
c < L 

is the effective modulus in peel, and the other quantities have been defined 

above. This constraint does not include any thermal effect terms because it 

derives from a purely plastic analysis and need not distinguish between how 

much of the adhesive strain was thermally induced and how much was developed 

by the applied mechanical loads. Were one to include the elastic stresses also, 

one could extend this thickness ever so slightly for very low thermal mismatches 

and/or very low temperature differentials. 


EXAMPLE 3: 

While this example should be regarded as only approximate because it deals with 
a very brittle resin system (and the theory above is precise only for reasonably 
ductile systems), it may be of interest to examine the permissible size of 
boron-epoxy infiltrated rods in aluminum extrusions. One uses equation (77) 
with an aluminum area equal to that of the boron-epoxy as a typical example. 

The resin has properties of the order of n = 0.005 inch, = 9000 psi , 

(y + Y ) = o.i, so that, if the hole in the extrusion be of diameter d and the 
temperature differential be Loo °F (from 350 °F to -50 °F), 

/ 4 4 \ 

2 x 9000 x Trd x 0.005 x 0.11 + I = [ (13.0- 2.7) x 10 -6 x LOO] 2 

\30 x 10 6 x ird 2 10 x 10 6 x Trd 2 / 
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whence d = 0.283 inch. 


It is evident that the 0.25 inch diameter selected by AVCO is close to optimum 
for the 350 °F cure-in-si tu resin system which has been used to date. To 
increase this diameter significantly requires either a decrease in the cure 
temperature or a layer of ductile adhesive between the boron-epoxy and the 
aluminum. AVCO is currently increasing the diameter of the reinforcing rods by 
use of a room-temperature-curing resin system which is subsequently post-cured 
to improve the strength and environmental resistance. This scheme is effective 
for two reasons. First, the thermal stresses are less severe in going from 
room temperature to either -67 °F or to the typically 250 °F post-curing temp- 
erature than is the case in going from 350 °F to -67 °F. Second, in the room- 
temperature-curing system, the stress-free state is at room temperature rather 
than at 350 °F. Therefore the residual stresses induced in the aluminum are 
far less severe and the fatigue life of the structure is improved. The major 
thermal stresses with this scheme occur at 250 °F and, as far as the bond is 
concerned, these are not subject to creep because there is no driving force at 
room temperature. 

The analysis above concludes, in equation ( 77 ), with an estimate of the maximum 
possible size of composite reinforcement which can be bonded to the basic metal 
structure. If mechanical loads t x and/or t 2 are applied at the end of the com- 
posite, where the maximum induced shear transfer stresses are located, this size 
will necessarily be reduced in order to accomodate both components of load. The 
understanding of this phenomenon permits of simple design modifications to all- 
eviate this difficulty. This is illustrated in Figure 8 in which, in the upper 
configuration, both mechanical and thermal loads combine at the same location. 

In the lower configuration, on the other hand, the extension of the reinforce- 
ment beyond the last fastener and beyond the softening cutout serves to separ- 
ate the peaks of the two transfer stress components. As explained in Section 4 , 
the peak mechanically-induced stresses are located at the last attachment and 
decay towards zero at the end of the member. The peak thermally-induced stress- 
es necessarily occur at the extremity of the reinforcement. Therefore, the 
stress-concentration relief shown in the lower part of Figure 8 enables the 
resin matrix to react only one peak load condition at any one location, permit- 
ting the application of higher mechanical stresses to the reinforcement. 
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3.2 Feathered-Edge Doublers 


For applications in which thick reinforcement sections are required, the abrupt 
cut-off at the end of the reinforcement is intolerable and the appropriate pro- 
cedure is to feather the end, as illustrated in Figure 9. In this event, quite 
a different mathematical solution governs and peel stresses (or interlaminar 
tension stresses in the composite reinforcement) cease to be a problem. The 
design of this latter class of joints is dominated by two or three prime con- 
siderations. First, the net extensional stiffness of the metal plus filament- 
ary composite parts should be constant throughout the joint, if at all possible, 
to promote maximum joint efficiency. Second, the net strengths at each end of 
the joint should not be so adversely unbalanced as to leave a long piece of 
reinforcement unable to accept any load because the end attachment area is too 
weak to load up the middle. Third, because the coefficients of thermal expan- 
sion of metals are distinctly higher than those of filamentary composites and 
because the stress-free temperature (just less than the cure temperature) is 
almost invariably higher than the operating temperature, a compressive load is 
usually more severe on the bond at each end than is a tensile load. 


The notation and sign convention for this type of joint are defined in Figure 
9. Again, deflections perpendicular to the axis of the filamentary reinforce- 
ment are excluded. Using classical mechanics of continuous structures, the 
equilibrium of the differential elements is given by 


- VT =0 
P 


+ WT =0 
P 


while the stress-strain relationships for the materials prescribe that 


cLx (EA) 


+ a AT 
c 
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+ a AT 
m 


( 83 ) 


in which 


d6 T 

m _ m 

dx (EA) 

m 

(EA) = E A (x/£) 
c c c 


(84) 


and 


(EA) = (E A ) j - [(E A )j - (E A ) 2 ](x/£) 
m mm 1 mm 1 m m z 

= (E A )i - A(E A ) (x/£ ) 
mm 1 mm 


(85) 


while displacement continuity requires that, for the adhesive. 


Y = 


(<5 - <5 ) / n 

c m 


( 86 ) 


The design problem may be expressed in its simplest form as the determination 
of the elngth, i, of scarf necessary for a uniform bond stress x to transfer 

Jr 

the load. Usually a compressive load is more critical than a tensile load 
because the thermally-induced residual bond stresses tend to relieve applied 
tensile loads and to aggravate compressive loads. In the process of solving 
the equations, it is necessary to check on the maximum adhesive shear strain 
induced since this serves as the failure criterion. The adhesive shear strain 
is assumed to be zero at the inboard end of the scarf. Even if it is not pre- 
cisely so there, it is zero close nearby. The numerically greatest adhesive 
shear strain will develop at the outer (feathered) end of the composite, regard- 
less of whether the applied load be tensile or compressive. Thus, the solution 
proceeds as 


dy 1 1 X w£ 

= -\(a, - a )AT + — 


(p - XpWX ) 


m 


dx nf w “ EA (E a ), - [(E A )y - (E A ) 2 ](x/fc) 

' cc mm 1 mm a mm" 1 

in which p is the total (tensile) load applied, whence 


x p v£ 


P£ 


Y = const. + — {(a -a )ATx + 

D 1 ° “ E A 


x + 


c c 


A (E A ) 
m m 


£n[(E A )j - A(E A ) (x/£)] 
m m ^ mm 


x w£ 2 


[A(E A )] 2 L 
m m 


(E A )j - A(E A )(-) - (E A ) j£n[ (E A ) 2 - A(E A )(-)] 
mm- 1 mm \ / mm- 1 mm- 1 mm 


(87) 


( 88 ) 


The substitution of the conditions y = o at x = z and y = y g + y p at x = 0 
serves to evaluate the integration constant. A further equation follows from 


30 



the assumed zero differential displacement across the adhesive layer in the 
(uniform) central region of the overlap. From equations (82) and (83), with 


6 = 6 , 

m c 


and, since. 


(T ) 2 (T ) 2 

2 - — = (a - a )AT 

(E A ) 2 E A ° m 


(T ) 2 = tvI and (T ) 2 = P - t vl , 
C p m P 


- T W& 

P 


E A (E A ) 2 
c c mm 


+ (a - a )AT 
c m 


Equation (91) permits the load p to be eliminated from the failure criterion 
formed with equation (88) applied at each end of the scarf. This process leads 
to an expression for i which may in turn be substituted into equation (91) to 
evaluate the maximum load capacity of the composite-reinforced structure. The 
equation for i becomes 


Ai 2 + B£ + C = 0 , 


where 


A = T w 
p 


,E A (E A )j - (E A ) 2 
cc mm mm 


(E A ) 2 £n[(E A ) 2 /(E A )i] 
m m z m m mm 1 

(E A )i- (E A ) 2 
mm mm 


/ (E A ) 2 ^n[ (E A ) 2 /(E A Ji3 

B = (a - a )AT 1 + 

c m \ (E A )i - (E A ) 2 

v mm 1 m m 


C * ^(*=0) - Ax=0 ] “ ±r,(v e + V f0r the upper b0und ) _ (95) 

= ±dy for the lower bound ) 

p 

The apparent problem of evaluating fcn(o) for the case ( E m A m )2 = 0 does no ^ 
fact arise because 

Oftn(O) = fcn(0°) = Jtn(l) = 0 . (96) 

Strictly this limiting case of a pure scarf joint is not covered by the present 

analysis because the thermal stress picture can change to the extent that y f 0 
at x = a. 
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The reason for the ± sign in equation (95) is explained in Figure 10. The 
mathematical confirmation of the precise adhesive shear behavior depicted in 
Figure 10 can be established by the following reasoning. Equation (87) may be 
re-arranged by means of equation (89) to read 

d ^ d ( ) _ (97 
dx ’ * (E mVA 1 

It is clearly evident that dy/dx = o only at and beyond x = z, except for the 
particular positive value of p, for each joint, for which dy/dx = o at all 
values of x. At x = 0, 


dy 

dx 


lfp[ (E A )i - (E A ) 2 ] - t w«,(E a )j 
_ mm 1 m m p m m _ 

n (E A h(E A ) 2 

L m m m m 


(98) 



a ) AT + 
m 


E A 
c c 


(E A h J 

mm 


(99) 


and dy/dx | o will tend to be positive as the result of usual thermal mismatches 
and compressive loads, so that y 0 will be negative while, for large tensile 
loads P, dy / dx | o will tend to be negative and y 0 positive. 


Returning now to equations (82) through (85) it is evident that A = 0 for 
e A c = A ( e^a^ ) which promotes longer effective scarfs as the result of carefully 
matching extensional stiffnesses along the length of the joint. (It is estab- 
lished in Reference 3 that only the unbalanced varieties of scarf joints suffer 
from limitations in joint efficiency.) In other words, it is sound design 
practice to make 


(E A )i = (E A ) 2 + E A 
mm 1 mm* cc 


( 100 ) 


The evaluation of the optimum length i for a design governed by compressive 
loads is straightforward. The length i is determined from equations (60) 
through (63) by reversing the sign of (a - <* m )AT and using the positive sign 
for y in equation (95). Actually, all the other quantities change sign except 
for the thermal terms, but it is simpler to compensate in this manner. The 
load capacity then follows from equation (91), again with the sign of the 
thermal stress terms reversed. The case of tensile loading is a little more 
tricky. The blind use of equations (82) through (85) with positive values of P 
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could easily result in a joint design which will hold its load only if applied 
gradually during cooling down in the autoclave after bonding. Once the load 
were released, the thermally-induced stresses could cause failure of the 
structure. In other words, for tensile loads it is necessary to check also the 
no-load condition to ensure that it is not more severe on the adhesive than is 
the load condition itself. For tensile loads, the appropriate sign in equation 
(85) follows from equation (99). If dy/dx|o is positive, yo is negative, and 
vice versa. The shear stress yg is usually negative for zero mechanical load. 

In the design of composite-reinforced structures it is necessary to guard 
against discontinuity problems in the transverse direction also. It is point- 
less to carefully taper the ends in the primary load direction and leave a 
severe thickness discontinuity in the orthogonal direction. An interesting 
discussion of this problem and its solution is given in Reference 7. 

The analyses above do not lend themselves to ready graphical representation 
because of the large number of variables involved. It is appropriate to demon- 
strate the conclusions and the use of the analysis by simple numerical examples. 
Consider, in turn, balanced and unbalanced joints to show the adverse effects 
of poor detailing and thin and thick sections to demonstrate how thick sections 
tend to exceed the predicted bond load capacity. Since the boron-epoxy to 
aluminum combination is receiving attention for selectively-reinforced struct- 
ures the examples discussed will be so oriented. 

EXAMPLE 4: 

For a bonded composite-reinforced metal structure of the type illustrated in 

Fiqure 10, let t = 0.100 in., (t )2 = O .060 in., (t )j = 0.360 in., w = 1.0 
3 ’ c m m 

in., E = 10X10 6 psi , E = 30X10 6 psi , a = 13.0X10 -6 /°F , a = 2.7*1CT 6 /°F, 
m c m. c 

at = - 50 - 350 = -Loo °F , Tp = 6000 psi, n = 0.005 in., and (y e + y p ) = 2 . 0 . 
Then, from equations (92) through (95), 

A = 0, B - (2.7 - 13.0)xl0 -6 x(-400)[l + ij£n(|) ] = O.OO 26 L 5 , C = 0.005*2.0 = 0.01 
whence, for compressive loading, 1 = -c/(-b) = 3.780 in.. From equation (91), 


33 



P = 10xl0 6 x0.06 


6000x3. 78OX 


- — )- 

O 6 x0.06/ 


(2.7-13.0)xl0 -6 x(-l+00) 


^30x10 6 x 0.1 10xl0 6 x0.06^ 

= 27,218 - 2,1+72 = 2UJU7 lb/inch width in compression and p = 27,218 + 2,1+72 
= 29,691 lb/inch width in tension. 

These loads correspond to composite stresses of 226,800 psi in both load senses, 


As the maximum allowable adhesive shear strain is reduced, the theory above 
predicts a progressive decrease in joint load capacity. The theory is conserv- 
ative with respect to brittle adhesives because of the somewhat surprising 
behavior of perfectly elastic scarf joints subject to thermal mismatch between 
the adherends. It is established in Reference 3 that, for both very short and 
very long overlaps (but not those in between) adherend thermal mismatch does 
not prevent the design of a scarf joint to transfer a required load between 
specified adherends. The dominant adverse effect is any stiffness mismatch 
which may exist between the adherends at each end of the joint. In the present 
context this scarf joint behavior translates into a statement that, provided 
the total extensional stiffness remains uniform throughout the joint region, a 
sufficiently small scarf angle can always be found to enable any adhesive, no 
matter how weak or brittle, to effect an adequate load transfer. The prime 
drawbacks to reliance on this idea are that, in the first place, the scarf may 
need to be so long that there is no room in the structure for a constant section 
between the joints at each end and that, secondly, the elastic analysis to 
determine such a scarf length is necessarily extremely difficult, being even 
more complicated than for a scarf joint (Reference 3). 


Example 5: 


Suppose that all data are the same as for Example 4 with the single exception 
of the presence of adherend stiffness imbalance with (t )j = 0.160 in., so that 
there is no bend in the reinforcement. The metal end of the joint is not built 
up adequately to carry the load which the filamentary composite reinforcement 
could otherwise develop. 

From equations(92) through(95) now 
/l 1 


A = 6000 


30xl0 6 x0.1 10xl0 6 x0. 


— j 1 + 0.6fcn(jy|) 


= - 0.00161+6 
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B = (2.7 - 13.0) xl0 _6 x ( -400 ) x 1 + 0.6£n(:p|“f = 0.001695 


C = 0.005 2.0 = 0.01 

Then, for compression 

0.001695 - /( 0.001695) 2 + 4(0.00l646)(0.0l) 

SL = — 

-2 X (0.001646) 


= 2.003 in. . 


The associated load capacity then follows from equation (59) as 


P = 10xl0 6 x0.06x 6000x2. 003( + ) - (2.7- 13.0)xlO 6 x(_400) , 

\30xl0 6 x0.1 10xl0 6 x0.06/ 

= 14,421 - 2,472 = 11,950 lb/in. for compression, 

= 14,421 + 2,472 = 16,894 lb/in. for tension. 

These values are only about half those predicted for the stiffness-balanced 
joint examined in Example 4. This comparison serves to emphasize the 
importance of careful design detailing to ensure maximum joint efficiency by 
minimizing any imbalances which'may be necessary. 


EXAMPLE 6: 

Consider a joint of the same materials as in Example 4 but with the adherends 
four times as thick. The plastic scarf length remains unaltered at 3.780 
inches and the composite stress consequently drops to 56,700 psi . The load 
capacity increases to 

P = 10xl0 6 x0.24 6000x3. 78 o( + ] - (2.7- 13.0)xl0 _6 x(-400) , 

\30xl0 6 x0.4 10x10^x0.24/ 

= 68,o4o - 9,888 = 58,152 lb/in. for compression, 

= 68,o4o + 9,888 = 77,928 lb/in. for tension, 

but the adherend materials are being used much less efficiently. (Actually, 
as discussed in Example 4, greater reliance should be placed upon the elastic 
capacity of the bond in a scarf joint, so this inefficiency prediction for 
thick sections is somewhat conservative.) 
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4. BONDED METAL REINFORCEMENT AROUND BOLT HOLES IN COMPOSITES 


The application of bolted joints in advanced composite structures depends upon 
several factors. On an individual basis, load transfer by bolts in plain holes 
in isotropic pattern filamentary composites (0 o A5°/90 o /-L5°) can be quite 
efficient. On a specific weight basis there is a strong case for not reinforc- 
ing the laminate with metal but for building up the laminate thickness as nec- 
essary instead. However, while this usually gives the lightest (and frequently 
the least expensive) bolted joint, it may be excessively bulky. There remains 
a case for analyzing bonded metal reinforcements around the bolt hole(s) for 
applications in which space does not permit a thicker laminate to be used. 

There are two possible locations for such doublers: within the laminate or 
externally. The two schemes are depicted in Figure 11. The internal doubler 
has twice the effective bond area of the external doubler but, in practice, is 
found to be associated with joggled fibers and poor quality laminates. External 
doublers are generally to be preferred because of the higher quality laminates 
and superior fit over the area to be bonded. The feathered edge shown in Figure 
11 is reserved for thick doublers. The thickness beyond which uniform doublers 
cannot be employed effectively is determined by peel -stress considerations 
explained in Reference 1 [see also Equation (79) here]. 


The analysis of external bonded doublers begins with a demonstration that, for 
uniform thickness doublers, the optimum location of the bolt hole is in the 
middle of the doubler. The nomenclature and mathematical model are defined in 
Figure 12. Equilibrium of differential elements in each of the four areas 
requires that, for one side of the joint. 


dT ! 

+ T = 0 , 

dT 2 

T = 0 5 

dx P 

dx 

P 

dT 3 

+ T = 0 , 

and 

dT 4 

- T = 0 

dx P 


dx P 


The thermo-elastic material equations are 


( 101 ) 
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46 k \ 


dx Vi 


+ a, AT 
k 


(k = 1,2, 3 , it) , 


( 102 ) 


while the adhesive shear strains are given by 

6 2 - <Si 


Yl2 ~ 


and y 34 = 


6 4 - 63 


(103) 


n n 

The boundary conditions apply largely to dy/dx. For the assumed fully-plastic 
adhesive stress state, with reference to Figure 12, 

k = e/l . (104) 

The critical strain developed in the adhesive occurs at the location x = £, 
while the adhesive is stress-free at x = 0. The solution proceeds along the 
path 


and 


dyi2 

1 

"2 T 2 

dx 

n 

E.t. 
" 1 1 

<1 2 Y12 

2 

1 

dx 2 

n 

E.t, 
1 1 

d Y 34 

1 

2 T4 

dx 

n 

E.t, 
1 1 

a 2 Y 34 

2 

1 

dx^ 

n 1 

E.t. 
1 1 


+ (a. - a ) AT 
1 o 


T = 

P 


+ ( a . - a ) AT 
1 o 


(105) 

(106) 
(107) 


T = 

P 


(108) 


The solutions of these differential equations are 


A 2 t 


Y12 ~ A12 + B 12 x + 


P 2 


"X 


(109) 


2 G 


and 


A T p 

Y34 = a 34 + B 34 ( x - e) + (x - e) 2 

2 G 


( 110 ) 


The application of the boundary conditions to equations (105), (107), (109) and 
(110) needs the evaluation of the adherend stresses due to thermal mismatch. 
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It is assumed that the overlap is sufficient to induce the maximum stresses, 
given by equation (70). 


T = - T. = (a. - a )AT / (■ 
O 1 1 o k 


2 + 2 


E t E.t 

o o i i 


■) 


at x = e 


(m) 


The adhesive shear strain at x = o will, by inspection, be zero for the maximum 
load case, being usually negative due to thermal stresses prior to application 
of load. The various boundary conditions are 

a 12 = o (from y 12 at x = o), (112) 

Bj 2 = (a. - a )AT/n (from dy^/^x at x = 0 ) , (113) 


b 12 + 


1 


(a . - a )AT 
P 1 o 

— e = 


n 


1 + 


(J— . + — ) ♦ \ 

\E t E.t. / \E t E.t./ J M'E.t. E^t J 

OC 11 OO 11 


E t 

11 o o 


(from dy 12 /dx at x = e). 


(114) 


A34 = + -^12 e + (A 2 t^/ 2 G)< 


(from Y12 = Y34 at x = e)> (115) 


B 


( a . - a ) AT 
i o 


34 “ 


1 + 


2 2 

+ 


\E t E.t./ YE t E.tJ 

OO 11 00 11 


Pe r + M . 

\E.t. E t / 


2P 


n£\E.t. E t 
11 oo 


nE t 
o o 


A 2 t 


B 34 + 


(a. - a )AT 2P 
— ( £ - e ) * 1 ° + 


(from dy 34 /dx at x = e), (116) 

(from dy34/dx at x = (117) 


n 


nE.t. 

i i 


and 


(y + y ) = A34 + B 34 ( £ - e ) + (X 2 t /2G)(fc-e ) 2 (from y 34 at X = t) . (118) 

e p P 

Elimination of A 34 and B34 by means of equations (111) through (117) leads to 
the following expression for the load carried: 


n(y + Y-) - (a. - a )AT£ 

p = T £ = — — i 

p U 2 \ n - 2e 


2 E.t. 
1 1 


)(vJ 


(119) 


The load is seen to be greatest when 

e = 1/2 . (120) 

Therefore, for efficient design, the bolt (row) should be located in the middle 
of the doubler. The adhesive shear strain at the middle of the doubler then 
follows from equations (83) and (88) as 


39 



Y - A34 - 


(a^ - a Q )AT/£ 


X^T p /£\ 2 


( Y e + V 


8 \E t 
o o 


( 121 ) 


which is precisely half that at the critical end (x = £) in the event that the 

doublers are balanced with respect to the basic adherend. That is E t = E.t.. 

0011 

While the stiffness ratio (E.t.)/(E t ) has no effect on the total load trans- 

11 OO 

ferred, it can be shown that the adhesive shear strain gradients are minimized 

when E.t. = E t . In this event, and with the bolt in the middle of the doub- 

1100 

ler, dy 12 /dx equals dY 3 i t /dx at x = i/2. Equation (119) seems to suggest that 

p decreases as 1 increases. Actually, however, 1 is not really a variable 

since the length of overlap is determined by the conditions that, at p , 

111 

Y 34 = (y 6 + y p ) at x = t and y 12 = 0 at x = 0 . Re-arrangement of equation (119) 
leads to the following expression for i 


g— j— i 2 + (a i - a Q )AT i - n(y e + y ) = 0 , 

i i * 


( 122 ) 


whence 


/ „ ^V^Ye + Y P^ 

-(a. - a ) AT + J[(a. - a )AT] 2 + — 

1 O W 1 O E . t . 

■ 11 


and, for each side of the joint, 


( 2 t /E.t. ; 
p 1 1 


(123) 


P = T £ = 
ult p 


/ ^ T Tj n (Yo + "O 

1 1 

2 /E. t . 

1 1 


- (a. - a )AT 
1 o 


(124) 


Strictly, this same load p given by equation (124) applies for any 

111 u lmate 

overlap > £ given by equation (123). The reason is the negligible load trans- 
ferred by the elastic adhesive in the additional parts of the overlap. In 
determining £ from equation (123), care must be taken to ensure that the maximum 
positive and negative loads are inserted in turn. If no reverse load is ever to 
be applied, the appropriate other limit is zero load to protect against a 
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possible design which would hold together as long as the load was maintained, 
but split apart as soon as it was unloaded. 

EXAMPLE 7: 

To illustrate the application of the formulas above, consider a graphite-epoxy 
composite laminate 1.0 inch wide and 16 plies thick of the (o 0 / 45 < 7 o 0 /- 45 °) 
pseudo-isotropic pattern. This is to be loaded up through bonded titanium 
doublers and a single bolt in double-shear. The composite adherend strength 
outside the joint area is 103,000 x l.o x 16 x 0.0055 = 9,064 lb. A suitable 
bolt to carry this load in double shear is a 5/16 inch diameter AN or NAS bolt 
heat-treated to 125 ksi. The net section strength of the laminate through the 
bolt hole is adequate for the half-load carried there. There is no need to 
consider the use of a smaller bolt of higher heat-treatment. Some of the bolt 
load will be transferred to the composite by direct bearing on the bolt but, 
because the modulus of the laminate is only about half that of the titanium, 
it is only slightly conservative to assume that all the bolt bearing takes 
place on the titanium. With a yield bearing stress of 216 ksi at e/d = 2 for 
annealed 6 a£-4v titanium, each doubler needs a minimum thickness of 0.067 inch 
for the 5/16 inch bolt. Then, with a ductile adhesive, for which x p = 6000 
psi , n(y + Y p ) = 0.010 inch, ou = 0, a Q = 5.8 x lO- 6 /°F, E_ [ = 11. 9 x 10 6 psi , 

t = 0.088 inch, at = -50 -350 = - 4 oo °F, the maximum load that can be carried 

1 

is evaluated by means of equation (124) as P ^ = 6,804 lb which, when doubled 
to account for the duplicated bond area, is adequate for the load required. 

The design minimum overlap is evaluated by means of equation (123) as l - 1.13 
inch and is 50 per cent more than the theoretical minimum for a perfect bond. 
This, then is a satisfactory design with a margin for environmental deterior- 
ation. A suitable actual overlap is 1.375 inch to allow adequate tolerances. 
Had the computed margin of potential bond strength to actual laminate strength 
(outside the joint) been less than 50 per cent, or even negative, it would have 
been necessary to increase the extensional stiffness of the composite locally 
by adding o° plies. It should be noted that the bond strength increase goes up 
slightly less rapidly than /Et. Also, it can be seen from equation (124) that, 
for thick laminates and brittle adhesives, there is a distinct possibility that 
the doublers will split off. In this case, checking equation (83) of Reference 
1 for peel stresses shows them not to be a problem. Assuming one bond-line 
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thickness of adherend affected by the peel stresses. 


/8000\ 4 l6 X10 6 X 0.005/ 1 1 1 \ 

t = ( } 1 + + 1 = 0.2^5 inch 

°max \6000/ 3 x(l-0.09) '0.5 X 10 6 1.7 x 10 6 l6xio 6 / 

while the actual value of t is 0.067 inch. 

o 
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5. IN-PLANE SHEAR TRANSFER THROUGH BONDED JOINTS 


The direct transfer of in-plane shear loads from one member to another is 
covered by the analyses in References 1 and 3. In essence, the governing 
equations are the same as for direct loads except that the adherend stiffnesses 
Et are replaced by the appropriate Gt. The thermal stress picture is more 
complicated, being a negligible problem except at the two diagonally opposite 
corners of the bond area because, elsewhere the thermally-induced adhesive 
shear strains are orthogonal with respect to those caused by mechanical loads. 
This situation is depicted in Figure 13. An effective simple technique for 
dealing with thermal stresses under in-plane shear loading is to use a 
modified form of equation (77) to compute the residual bond strain as 


Y 


r 


[ (<*2 



(125) 


and subtracting this from the total adhesive shear strain otherwise available 
for reacting the mechanical loads. 


The analysis of a single load path for in-plane shear is straightforward, but 
that for multiple load paths, as at the intersection of two cells of a multi- 
cell torsion box is more complicated. An illustrative example of such problems, 
below, serves to describe the technique to be employed in such cases. Thermal 
stress considerations have been omitted because, in practice, such structures 
are usually made of a single material throughout. 


The mathematical model of the structure to be analyzed is defined in Figure 14. 
The governing equations, for the sign convention adopted, are as follows: 

as i 

+ T = 0 , 

dx 

as 2 

- T = 0 , 

dx 

as 3 

+ T = 0 , 

ay 


d<S 


1 


>1 


dx 

d6 2 

dx 

as 3 
ay 


Git i 

S 2 

G 2 t 2 

s 3 

G 3t3 


(126) 

(127) 


(128) 
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dS 4 

d 6 4 

s 4 


- 2x = 0 , 


” 9 

(129) 

dy 

dy 

G 4 t 4 


dS 5 

d 6 5 

s 5 



+ T = 0 


dy 

dy 

G 5 t 5 

dS 6 

“ T = 0 , 

d6 6 

s 6 

dxi 

dxi 

G 6 t 6 

dS 7 

+ T = 0 , 

d 6 7 

S 7 

dxi 

dx^ 

G 7 t 7 


The solutions of these differential equations are 

(^12) 2 T 


Yl2 ~ A 12 + B 12 X + 


Y34 _ a 34 + B 34 y + 


Y45 “ a 45 + B 45 y + 


P 2 

—DC 


2 G 

C A 3J+ ) 2 t 


2 G 


(X4 5 ) 2 i 


2 G 


P 2 

-y 


-y 2 


(130) 

(131) 

(132) 


6 2 - 

4 1 


«4 

- 5 3 

«4 - 

«5 

6 e - 67 

(133) 

Y 12 ~ 

n 


9 

Y 34 ~ 

> Y 4 5 - 

n 

n 

9 

Y 67 = 

n 

It follows 

that 







dYl2 

h 

f S2 

s > \ 

d 2 Yl2 T p 

( 1 

1 1 

| (^12 ) 2 

(134) 

dx 

d 

^G 2 t 2 

Gitj/ 

dx^ n 

Wi 

+ 

G 2 t 2 ^ 

, = ’ 

dY 34 

h 

i s 4 

S 3 V 

d 2 Y 34 T pj 

( 1 

2 > 

1 

\ (^34 ) 2 

(135) 

dy 


^4^4 

G 3 t 3 / 

r- 

1 

04 

& 

#1 

(c 3 t 3 

G 4 t 4 ^ 

,= 0 T p ’ 

dy 4 5 

1 

_| 

( ^4 

S 5 v 

d 2 Y 45 Tpi 

< 2 

1 > 

I U-s) 2 

(136) 

dy 

1 

n 

'G 4 t 4 

G 5 t 5 

r- 

II 

Csl 

>> 

TJ 

n 

^G 4 t 4 

Gst 5 ^ 

•' G P ’ 

dYe 7 

l 

_j 

( s 6 

S 7 V 

d 2 Ye 7 

( 1 

1 ) 
— L 

1- <>67>2 t 

(137) 

dx 1 

! 

n 

^6^6 

G 7 t 7 

» O 

dx 2 2 n 

^6^6 

G 7 *b 7 ^ 

0 p ' 


(138) 


(139) 


(140) 
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Y67 = A 67 + B 67 x 1 + 


(141) 


2G 

in which x and x T differ only in origin. Most of the boundary conditions are 
applied by means of the first of each of the equation sets (134) through (137). 
Thus 


Si 

Bj 2 = - 5 at x = o , 


(*12 ) 2x t 


B 12 + 


'12 


T p fc 12 (Si - T p A 12 ) 


l_G 2 t 2 


Git i 


at x = 


(142) 

(143) 


B 


34 


b 34 + 


V 12 

nG 3 t 3 

(X 34 ) 2 T 


-34 - 


B 4 5 = - 


B 4 5 


T p £ 45 

nG 5 t 5 

P. 


(A45) ? ‘ t 


^45 - 


t d £ 34 


nG 3 t 3 

S 4 

r|G4t4 

t P ^67 

nG5t 5 

s 4 

nG 4 t 4 


at y = 0 , 


at y = £ 3 4 , 


at y = 0 , 


at y = £ 4 5 , 



(*67 ) 2t t 


b 6 7 + 


G 7 t7 


"*-67 = ~ 


G 6^6 J 


at X! = O , (x = & 12 ) 


at X 3 = &67 , (x - l \2 + £ 67 ) 


nG 2 t7 


(144) 

(145) 
046) 

(147) 

(148) 

(149) 

(150) 


In these equations, t 12 , t 34 , £ 45 and £ 6 7 re f er to the effective ful ly-plastic 
zones out of the total overlaps, which may be greater but do not contribute to 
extra load transfer. Since the angles tying the skin and web together may also 
serve as concentrated spar caps, none of the frequently possible reductions in 
variables is effected here. Nevertheless, since the outer legs of the angles 
are stress free, it follows quite generally that 


45 



*12 = *34 


and 


*■45 = *67 • 


(151) 


Gross equilibrium also requires that 

(*•12 + *67^ = ( s i ~ ^7 ) /Tp = S 4/ T p 


while, also. 


Now, from equations 


(*34 + £ 45 ) - S 4/ T p 

(144) through (147), 


(*•4 5 



G4t4 J 


0 


(152) 

(153) 


(154) 


so that, quite generally, 


*12 = *-34 = *45 = *67 - ^effective 


S 1 ~ S 7 


2x 


2t 


(155) 


P P 

In other words, each bond transfers precisely the same load as the other three. 


One of these four extents of overlaps becomes critical first, depending on the 
relative adherend stiffnesses. The maximum possible extent of plastic adhesive 
zone occurs when the shear strain is zero in the middle (strictly y g to 
identify between upper and lower bounds) and at its maximum value (y g + y ) at 
each extreme of the overlap. However, this is likely to occur only for the 
central web if G 3 t 3 = G 5 t 5 = G 4 t 4 /2. In general, one side of the overlap will 
be more critical than the other and the extent of effective overlap less. The 
derivation of the critical overlap follows from the scheme outlined in Figure 
15. Quite generally, the effective overlap is 


where 


and, here. 


£ = 




dy a /to \ Ft 
dy^/d x'V A 


dy 
* a 

< 

dy e 

dx 


dx 


A = ( W2G 


and y 6 = (Y e + Yp) 


(156) 


(157) 


(158) 
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Note that the minus sign in equation (156) automatically takes care of the 
possibilities that the gradients have the same or opposite signs. 


For the load transfer from element 1 to element 2, 


dy^/ dx 

dy /dx 
P 


Veff 

G 2 t 2 


(S.-Veff) 


Git 


1^1 




( 159 ) 


or, if this ratio is numerically greater than unity, its inverse. Likewise, 
between elements 3 and 4, 


dy /dy 
g 

dy^/dy 

and, from 4 to 5, 


G 3 t 3 \/ S 4 


V( W Veff 


(or its inverse) 


(160) 


dy^/dy 

dy^/dy 

while, from 6 to 7, 


(G 5 t5\/ S 4 

W 4 t 


(or its inverse) 


(161) 


dy /dx 
a 

dy 3/ dx 



(162) 


or its inverse, depending upon magnitude. In applying equation (156), it must 
be remembered that A is different for each zone because of the different values 
of (A ) 2 . Thus, the identification of the minimum effective overlap reduces 

g$ 

to the problem of evaluating the least value of 
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Once the minimum value of a J .. has been established in terms of the 

effective 

adherend properties and load ratios and the adhesive properties, the actual 
adhesive and adherend loads at failure are proportional to 


Su = 2 t 


"effective minim-urn 


(168) 


It can be shown easily that, to maximize the total load capacity of the bond, 
it is necessary that the overlapping adherends be as stiff as practical, even 
if this means a thicker member reduced in thickness away from the joint. This 
problem is too complex to lead to such specific recommendations that the adher- 
end stiffnesses should be balanced, except for the 3-4-5 elements for which, 
if the bond limits the total load capacity, maximum efficiency requires that 


^ 3^3 = ^565 = G 4 W2 


(169) 


The actual use of this analysis is illustrated in the following example. 


The analysis above is based on a hypothetical perfectly-plastic adhesive. 
Effectively, for real adhesives also, the plastic portion of the stress-strain 
curve contributes to most of the load transfer and it is only slightly conserv- 
ative to neglect the elastic load transfer in comparison. On the other hand. 
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it is important, in sizing the overlaps, to add to the computed plastic zones 
a sufficient distance for the elastic stress to build up from zero to the 
plastic value. Typically this distance is of the order of 3/ A so that, if 
dy/dx has the same sign at each end of an overlap, a distance equal to 6 /A 
should be added to the sum of the two plastic zones. 

EXAMPLE 8: 


Consider the structure and load ratios shown in Figure 16. Use the notation 
of Figure 14, with adhesive properties = 6ooo psi , y g = o.i, y^ = 2.0, and 
n = 0.005 inch. Assume a 2024-T3 aluminum structure for the skins, web and 
angles, so that g = i+.Oxio 6 psi and F oii = Lo ksi. Then the maximum effective 
overlap is the least of 


SU 


'12 


V 


2 x 0 . 0 Q 5 
6000 


x 2.1 | 


( 1 - 0 . 5 ) 




>V 

*-34 = ^45 = 3 . 7^2 x 


1 + , 1 

k x 10 b x 0.125 L x 106 x 


O. 25 )' 


= 1.350 inch 


»v 


1 + 2 X 


0.063 

0.125 


2 1 i..: 

0.125 0.063 


1 _ 


= 1.331 inch 


1 - 


l + 


'67 


= 3 . 7 ^ 2 * 


*2 0 . 063 / 


= 0.598 inch 


V (0.25 + 0.063 ) 

Thus the skin-to-angle bond governs the load capacity and (* effective ) minimum 
is equal to 0.598 inch. In turn, it follows that 

S 4 = 2 x 0.598 x 6000 = 7,176 lb/inch , 


Sj = 2 x S 4 = lL ,352 Ib/inch , 

S 7 = S 4 = 7,176 lb/inch . 

The corresponding shear stresses in the skin and web are at a uniform value of 


t = 7,176 x 16 = 11/4,816 psi 

which is beyond the allowable f , so the adherend strengths govern. The 
recommended value of a then becomes 
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n 67 = i4Qj g Q , ° ' 6ooo° 6 ' 3 * 1,5 + + °- 25 = °' 315 + °* 6T1 + °- 25 = 1 - 236 inch 


for the 6 to 7 overlap, with the 1.5 factor being the design margin of 
potential bond strength over adherend strength, 3/x the single elastic shear 
stress build-up zone, and the 0.25 the manufacturing tolerance. In the 1 to 2 
and 3 to 4 (or 4 to 5) zones, the recommended overlaps are, respectively, 

sl 12 = 0.315 + 6/A 12 + 0.25 = 0.315 + 1.732 + 0.25 = 2.297 inch , 

£ 34 = * 45 = o.3l5 + 6/A 34 + 0.25 = 0.315 + 1.225 + O .25 = 1.790 inch . 

Note that, even though the low-stress elastic troughs make up a large portion 
of the total overlaps, they contribute but little to the static joint strength. 
Their major important role is to attain a long service life by providing an 
area of inevitably small adhesive strain to resist creep and to serve as the 
plastic zone late in the life of the joint when the outer edges of the original 
plastic zones have deteriorated through environmental exposure. 
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6. PEEL STRESS RELIEF FOR DOUBLE-LAP JOINTS 
(TAPERED-LAP JOINTS) 


Reference 1 suggests certain simple design modifications to extend the thick- 
ness range over which double-lap joints can develop adequate efficiencies. 
Basically, the usual strength limit is due to excessive peel stresses which 
develop at the ends of the outer adherends at a load level which may be signif- 
icantly less than either the potential shear strength of the bond or the adher- 
end ultimate strength. Beyond a certain adherend thickness, for a given comb- 
ination of materials, the joint efficiency drops to intolerably low levels. 

The object of this section is to explain how to effect a major relief of the 
peel stress problem, thereby extending the economy in fabrication of uniform 
lap joints to a greater usable range of thicknesses. The design technique 
employed is that of tapering the ends of the outer adherends, as shown in 
Figure 17. Actually, the joint is slightly stronger if the thickness at the 
end is not quite reduced to zero. A tip approximately 0.010 inch thick is both 
easier to fabricate and handle without damage. This tip must be minimized to 
prevent the build up of peel stresses. The analysis below, for a completely 
feathered end of the overlap, is slightly conservative with regard to the ult- 
imate bond strength but is only slightly so and effectively eliminates one 
variable from the analysis/design process. 

A general analysis prepared for this problem proved to be so complicated as to 
be beyond algebraic manipulation without assistance from a digital computer. 
Therefore, what follows is a simpler procedure governing the approximate design 
of optimum proportions for static ultimate load only. Partial load conditions 
are not covered by this simplified design procedure but this is considered to 
be acceptable because such load conditions do not govern the geometry of the 
design. A slightly conservative analysis for the tapered-lap joints at any 
load level can be provided by ignoring the taper and using adhesive shear 
analysis for uniform lap joints. This particular analysis had been approached 
with a willingness to sacrifice some of the unrealized bond shear strength 
potential in order to develop the remainder. It transpired that it was quite 
simple to attain the entire potential bond shear strength while eliminating the 
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peel-stress problem completely. Indeed, by skillful design, it is possible to 
exceed by 24 per cent the shear strength potential of uniform double-lap joints 
by making each of the tapered outer adherends 33 per cent stiffer than one half 
of the uniform inner adherend. 


The mathematical model for the analysis is illustrated in Figure 17 along with 
the element forces and displacements. Thermal mismatch effects are excluded to 
elucidate the dominant influences. The small load transfer throughout the 
elastic trough is ignored and an analysis for a perfectly-plastic adhesive is 
conducted. At ultimate load the optimum design will maximize the sum of and 
s. the extents of the adhesive plastic zones in the uniform and tapered zones, 
respectively. The analyses of the zones A-B, B-C and C-D in Figure 17 must be 
performed separately and the solutions joined by matching boundary conditions. 
The location of C is defined to occur at the start of the taper on the outer 
adherends and the length sl is treated as an unknown to satisfy this condition. 
(This considerable simplification in analysis eliminates the matching of cond- 
itions at one more transition in behavior and is one reason why the analysis is 
inapplicable to partial load levels. Once designed for ultimate load, the 
start of the taper remains fixed but the width of the plastic adhesive zone 
varies with load level. A further reason is that the length would change 
with load direction in the presence of any adherend thermal mismatch.) 
Throughout the region A-B, the adhesive shear stress is constant at T p while 
the adhesive shear strain y varies, reducing to zero at B because the adherend 
strains are identical throughout the region B-C. (The same is true in the pre- 
sence of any adherend thermal mismatch and ensures that dy/dx = o at each end 
of the zone B-C). Because of the identical uniform adherend strain assumed in 
the zone B-C for all adherends, it is necessary that 


whence 


2 t p\ 


2t_4 
P u 


E t E.t. 

O O 11 


(170) 


E.t. I 

(171) 

E t l. 
oo t 


(Were precise allowance to be made for the elastic adhesive shear transfer in 
the region B-C, equation (171) would still apply at that single location for 
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which dy/dx = o ). Throughout the zone A-B, 


+ 2t =0 , 

P 


- 2x =0 , 


where T represents the force per unit width in both of the outer adherends. 
The adherend displacements are given by 


dS T 
o _ o 

dx E t 
o o 


dx E . t . 

l l 


while the adhesive shear strain is defined by 


Y = (<$£ - <5 Q ) / 0 • 


It follows that 


T 0 = P - 2 V 9 


T. = 2t x 
i P 


where p is the applied load, so that 

dy 1/ T. T 


dy 1/ T. T 

— = -f _ -2- 

dx n'E.t. E t 
11 o o 


whence 


2 \/x 


Y = Y 0 - 


x / 2 


r)E t n YE.t. E t / ' 2 

O O 11 o o 


y being the adhesive shear strain at x = o. Since a is defined to be such 
’o 3 u 

that y = o (or Y = Y e for a lower-bound solution), 


Y 4t U + On- T 
' o P u t u , 


ny E t 
e o o 


T / 2 2 \ 

-4— * — 

ny VE.t. E t / 
e ii o o 


l 2 = 0 . 

U 


At maximum possible load capacity, y q = Y e + y p - In the zone C-D, equations 
(172) and (174) still hold, but equations (173) and (175) become 


d<$ T 
o o 


dx E t E t U/l. ) 
o o o o t 


d6 . T. 
l _ l 

dx E . t . 

i l 


T = 2t c , 
o 9 


T. = P - 2t p ; . 


Hence 
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dy 

ae 


dy 

dx 


2t £ . 2t C P 
P t + £2 _ 


E.t. E.t. 
11 11 


E t 
o o 


(181) 


Using the condition y = y at c = 0, this equation may be integrated to read 


Y = Y £ + - 


2t £ 

P t 


PC 2 t 

+ — 


E t 
o o 


E.t. 

l l 


E.t. 

l l 


£) 


(182) 


and, with the condition y = o at c = £ , this becomes 


1 

bTl.Z. 

p t t 


- E t 


Ut (£ + £, )£, 

P v u t' t 

E.t. 

l l 


2tX£ ' 
P t t 


E. t. 

l i 


= 0 


(183) 


For the most efficient design, y will reach (y e + y p ) at the same load level as 

that at which y does the same at the other end of the joint. This ensures the 

o 

maximumization of the sum (& +z ). The matching of y and y. requires that 

tl O JC 


bZZ b(l + Z )Z 2 £ £ 

1 X» U U 1 ^ 1 

E t E.t. E.t. 

O O 11 11 


U(£ + £j£ 

u t 


E t 
o o 


E t / U 


^ E . t 

l i 


E t 
o o 


u 


(184) 


The inclusion of the constraint (171) reduces the number of variables so that 



It is easily shown that the solution of this cubic is 


(185) 


£ E.t. 

— = -1^= 1.321+7 . (186) 

SL E t 
u o o 

That is, the net extensional stiffness of the outer adherends should exceed 
that of the inner adherend by 32.5 per cent in an optimally designed tapered- 
lap joint. It seems to be significant that current metal practice in bolted 
joints of similar configuration calls for each of the (tapered) outer adherends 
to be between 25 and 40 per cent thicker than half of the (uniform) inner 
adherend. 


Because of this stiffness mismatch, the adhesive plastic zone A-B extends 
further than would be the case if e t = E.t.. For uniform thickness outer 

O O 11 

adherends with E t > E.t , the end D would be critical, with reserve capacity 

O O 11 

at end A. Returning to equation (178), and incorporating equations (171) and 
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(186) for the optimized design, it can be shown that 


ny E.t. 
e i i 


■U ) 2i 

u 


/I + 1.3247 \ 
\2 x 1.7247/ 


(187) 


Hence 


£ = 1.0676 
u 


2( W 


(Ux p )/(riY e E i t i ) 


(188) 


and is 1.07 times as great as would be the case if E Q t Q = e^ for a given 
inner adherend to be bonded. The effect of optimizing the design for peel- 
stress relief is, consequently, to increase the potential shear strength to 


P 


optimum 


1.0676 X 



1 . 3247 
2 


■) p 


"balanced 


(189) 


which is 24.09 per cent better than the best that could have been achieved 
with uniformly thick adherends. 


It seems remarkable that this relatively inexpensive peel-stress relief modif- 
ication should be associated with such a significant increase in potential 
bond shear strength. This simple design/fabrication modification permits the 
more economical lap joints to be employed effectively throughout a significant 
range of thicknesses which would otherwise have required the use of more 
expensive scarf or stepped-! ap joints. The influence on joint strength of this 
concept is depicted in Figure 18 for HTS graphite-epoxy. The tapered-lap 
joints without the extra build-up on the outer adherends are found to have the 
same potential shear strength as for uniform outer adherends. The end A is 
more critical than the end D (Figure 17) because both the peel and shear stress 
concentrations are relieved by tapering. That tapering the outer adherends in 
accordance with Figure 17 has no effect on the shear transfer when E^ = E Q t o 
follows from equation (171). 
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7. CONCLUSION 


This report shows how to apply the simple concept of shear load transfer 
through definable ful ly-effective zones to adhesive-bonded joints of non- 
standard configuration which occur frequently in aerospace practice. The tech- 
niques have been illustrated through a variety of examples which are by no 
means the limit of applicability of the method. Actually, because the solution 
of the governing differential equations is far simpler for the all-plastic 
adhesive than for the perfectly-elastic adhesive, the use of upper and lower 
bound plastic solutions expands greatly the range of problems for which closed- 
form analytical solutions can be derived. 

The adoption of this plastic zone approach offers a distinct advantage over the 
older method of computing bond strength as the product of a uniform allowable 
shear stress and the total bond area. The fact that most of the load is trans- 
ferred through the plastic adhesive zones has been established by more precise 
elastic-plastic analyses (References 1 and 2). Likewise, it was shown that the 
extent of such plastic zones is independent of the total (long) overlap. 
Therefore, the use of the plastic zone analysis technique correctly relates the 
potential bond shear strength to the appropriate elastic and geometric para- 
meters instead of to the total overlap area, of which precise analyses have 
shown the maximum bond strength between given adherends to be independent. 
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SIGN CONVENTION 


LOADS ON ELEMENT 
OF REGION 2 


LOADS ON ELEMENT 
OF REGION 3 


FIGURE I. CO-ORDINATE SYSTEM AND DEFORMATIONS IN ADHESIVE-BONDED 
DOUBLERS 
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SKIN EFFICIENCY Tj 



FIGURE 2. SKIN EFFICIENCY FOR BONDED DOUBLERS WITH SIMPLY-SUPPORTED 
EDGE CONDITIONS 
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SKIN EFFICIENCY 7 ? 



FIGURE 3. SKIN EFFICIENCY FOR BONDED DOUBLERS WITH BUILT-IN EDGE 
CONDITIONS 
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FUSELAGE 


REFERENCE: 
AFFDL-TR-66-143 (1966) 




SKIN JOINT, STRINGERS CONTINUOUS. 
2 STRAP PLATES OF SKIN JOINT. 

3 DOUBLER PLATES, ADHESIVELY BONDED 
TO SKIN. ONE DOUBLER IS MADE 
AS A FINGER PLATE. 


CRACKS COVERED BY PRESENT ANALYSIS 
DESIGNATED BY c. NOTE DISTINCTION 
BETWEEN THESE CRACKS AND THOSE 
ASSOCIATED WITH 
MECHANICAL ATTACHMENTS a. CROSS-SECTION A - A 



0 



CROSS SECTION 
AT RIB 0 

FRONT 
SPAR 
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CO 


REAR ijfl 
SPAR > 7 




° 1 

(r} 

o o \ 

\ 
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c | 

FINGER ' 

SKIN 
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o o 

DOUBLE Ri 
o o / 


A 


FIGURE 4. FATIGUE CRACKS DEVELOPED DURING TESTING OF AIRCRAFT 
STIFFENED PANELS WITH BONDED DOUBLERS 
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s 


w 


OUTER EDGE 


PLASTIC 


ELASTIC 


U* 


ADHESIVE BEHAVIOUR 


INNER EDGE 
1 — d 


! L 


-SHEAR MODULUS G |E d x / | 
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RIVET t FOR INNERMOST ROW OF RIVETS 


RIVET FOR POSSIBLE ADDITIONAL ROW OF RIVETS 
TYPICAL CROSS-SECTION 
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Ax (REFERENCE) 


CD 


T « 




1 ♦ — — Ax 

E s t s ; 


6 s ♦ 


Ax (REFERENCE) 


DISPLACEMENTS AND FORCES 


FIGURE 5. CO-ORDINATE SYSTEM AND DEFORMATIONS IN BONDED DOUBLERS 
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JOINT GEOMETRY 




STRESSES 



DISPLACEMENTS AND ELEMENT FORCES 


FIGURE 7. NOTATION AND MATHEMATICAL MODEL FOR ANALYSIS OF 

TRANSFER STRESSES WITH SELECTIVELY-REINFORCED STRUCTURES 
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STRINGER RUNOUT CONFIGURATION TYPICAL OF METAL PRACTICE 



RESIN MATRIX SHEAR STRAIN 
DUE TO RESIDUAL THERMAL 
STRESSES 


RESIN MATRIX SHEAR STRAIN 
DUE TO MECHANICAL LOADS 


NOTE HOW PEAK SHEAR STRAINS COINCIDE ABOVE FOR EACH LOAD COMPONENT 




STRINGER RUNOUT CONFIGURATION SUITABLE FOR SELECTIVE REINFORCEMENT 

BY INFILTRATED EXTRUSIONS 



RESIN MATRIX SHEAR 
STRAIN DUE TO 
THERMAL EFFECTS 



RESIN MATRIX SHEAR 
y STRAIN DUE TO 
4 ) MECHANICAL LOADS 


NOTE HOW PEAK SHEAR STRAINS FOR THE TWO LOAD COMPONENTS 
ARE SITUATED AT TWO DIFFERENT LOCATIONS 


FIGURE 8. DESIGN DETAILING FOR SELECTIVE REINFORCEMENT OF METAL 
STRUCTURES BY BONDED COMPOSITES 
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BOND WIDTH w PERPENDICULAR TO PAGE 
JOINT GEOMETRY 


COMPOSITE 
“ STRESS 



METAL 

STRESS 




(EA) c = E c A c (x/£ ) 

(EA) = (E A ), - [(E A ). - (E A ) 0 ] (x/l) 
* 'm ' m m'l LX m m'l ' m m'2 1 ' ' 


ELEMENT FORCES 


FIGURE 9. NOTATION AND MATHEMATICAL MODEL FOR ANALYSIS OF 
BONDED DOUBLERS WITH FEATHERED EDGES 
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JOINT GEOMETRY 



FIGURE 10. ADHESIVE SHEAR STRAIN DISTRIBUTION FOR BONDED TAPERED 
DOUBLERS 
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INTERNALLY -BONDED DOUBLERS (IDEAL) 



VOIDS IN LAMINATE 



LONG CRACK 

BLUNT LIPS ON ENDS OF DOUBLERS 
(JOGGLED PLIES IN COMPOSITE) 


INTERNALLY -BONOEO DOUBLERS 
(REPRODUCED FROM ACTUAL PHOTOGRAPHS) 



EXTERNALLY- BONDED DOUBLERS (ACTUAL) 


FIGURE II. INTERNAL AND EXTERNAL BONDED METAL DOUBLERS 
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! ! i 

(NOTE THAT ENTIRE BEARING LOAD IS ASSUMED TO BE REACTED ON METAL 
DOUBLER(S) BECAUSE OF HIGHER MATERIAL STIFFNESS THAN FOR COMPOSITE). 

LOAD DISTRIBUTION 


FIGURE 12. NOTATION AND MATHEMATICAL MODEL FOR ANALYSIS OF BONDED 
METAL DOUBLERS AT BOLT HOLES IN COMPOSITE LAMINATES 
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BOND- LINE DISPLACEMENTS DUE TO THERMAL MISMATCH 



DIRECT LOADING IN-PLANE SHEAR LOADING 


(>' e + 7p) 2 = (7 r ). 


O'e + V 2 = l ( ’V )/| 2 


+ ^mechanical + 


y • , + / y } 

mechanical v/ r y ’ 


y , . , = MAXIMUM ADHESIVE SHEAR STRAIN 

mechanical |NDUCE[) By MECHAN1CAL LOADS 

y = RESIDUAL ADHESIVE SHEAR STRAIN DUE TO THERMAL 
r MISMATCH BETWEEN THE ADHERENDS 


FIGURE 13. ADHESIVE SHEAR STRAIN MAXIMA UNDER MULTIPLE LOADS 
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LARGE ARROWS INDICATE DIRECTION ASSUMED FOR LOAD TRANSFER 


SIGN CONVENTION 


FIGURE 14. NOTATION FOR ANALYSIS OF SHEAR FLOWS IN MULTI-CELL 
TORSION BOXES 
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FIGURE 16. DETAIL OF MULTI-CELL BONDED TORSION BOX 
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IDEALIZED ADHESIVE SHEAR STRESS DISTRIBUTION 


FIGURE 17. NOTATION AND MATHEMATICAL MODEL FOR ANALYSIS OF 
TAPERED-LAP BONDED JOINTS 


JOINT STRENGTH P (LB/IN.) MAXIMUM POSSIBLE JOINT EFFICIENCY T\ 


FIGURE 
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18. COMPARISON OF JOINT STRENGTHS FOR TAPERED AND UNIFORM 
DOUBLE-LAP JOINTS 
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TABLE OF MATERIAL PROPERTIES FOR FIGURE 18 


HIGH-STRENGTH GRAPHITE-EPOXY: 

(0 o /+45 o /90 o /-45 o ) s pattern: 

e£ = 8.0 x 10 6 psi, E^ = 1.7 x 10 6 psi, 

F^ u = 69 ksi , F^ u = 8 ksi , 

(in which the subscript N refers to properties in the thickness direction). 
DUCTILE ADHESIVE: 

t = 6 ksi , n = 0.005 in., y h ~ 20, 

p p e 

n (i +Y ) = 0.0102 in., e - 500 ksi, a - 10 ksi. 

u 2'e 'n C c 

^ c r max 

BRITTLE ADHESIVE: 

t = 9 ksi, n = 0.005 in., y/y a = 1-5, 

P P e 

n (i Ye + Y ) = 0.00042 in., E c = 1500 ksi , ct c =17 ksi. 





APPENDIX 


The following Fortran IV computer programs were developed in the preparation of 
Figures 2 and 3 for the efficiency of skin material when used in conjunction 
with bonded doublers around the edges on one side only. No data cards are 
required since the parameters are read in as part of the program. Illustrative 
sheets of output are included for each program. The first program applies for 
built-in edges at the outer periphery of the skin/doubler combination while the 
second one covers simply-supported edges. The doubler *./t parameter is defined 



and the bending stiffness parameter as 


bendk(m) = 
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PTMPMcr n\» m (38), S IGAVf 3« * 10 ) • FSIGMMIO)* °FNPK( 13) 
r sf t OVE °l AP A° p av 


C SET 


mm *■ 

0 

. 

OL ( 2 ) = 

0 

. 1 

ni ( 3) = 

0 

. 15 

OK 4) = 

0 

.2 

PL(5) = 

0 

.25 

HL ( 6 ) = 

0 

. 3 

nt_( 7 ) = 

0 

.35 

nt(8) = 

0 

. 4 

01(9! = 

0 

. 5 

OHIO! 

= 

0.6 

nuin 

= 

0.7 

M C 1?) 

= 

0.8 

PK 13) 

3 

0.9 

QK 14) 

= 

l. 

ni< 15 ) 

= 

1.5 

01. ( 16) 

- 

?. 

01(17) 

* 

2.5 

OLI1R! 

= 

3. 

nil 19 ) 

= 

3.5 

01(20) 

= 

4. 

OL ( ? 1 ) 

= 

5. 

OK 22) 

- 

6. 

0L< 23) 

= 

7. 

OK 24) 


8. 

OK 25) 

= 

9. 

11(26) 

- 

10. 

OL ( 27 ) 

= 

15. 

OK 28) 

= 

20. 

PL( 29) 

= 

25. 

OK 30) 

= 

30. 

PK 31 1 

= 

3 5. 

OK 32) 

= 

40. 

OL ( 33) 

= 

50. 

OL( 34) 

= 

70. 

0L( 35 ) 

= 

100. 

OK 36) 

= 

200. 

OK 37 ) 

= 

500. 

0L< 38 > 


1000 


PARAmcjfp ARRAY 


RFMOKf I > 
RENDK< 2 ) 
F^PNJPK ( 3 ) 
RFNDK ( 4 > 
RF^ni<{ 5 ) 

3F*nK< 6 ) 

8 ENDK ( 7) 

R ENnK ( 8 I 
R ENOK ( 9 ) 
RFNDK ( 10 ) = 


oa 

0. IS 
0.2 
0.3 
0.5 

0. 7 

1. 

1.5 

2 . 

3. 


REN^M 1 1 ) = 5. 

BENPK112) = 7. 

R ENOK (13) = 10. 
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Dinvn4 = T IPVT 4 V* 3 

C SFT TNI T T 41 ESTIMATE ON AVE a AGP A DHF&FNO STRESS 
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P = V3 / V4 

IF C (1.0001 .GT. P ) .AND. (0.9999 .LT. R) ) GO TO 20 
V3 = V4 
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30 CONTINUE 
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3 3X • 7HD0U R L F R , 18X, 4 Q HA VFR AGE ADHFRFNO STRESS / MAXIMUM AnuFREND 

4 STRESS / » 5 X , 3HL/ T /, 3X, 9H P A R AM F T F 0 , 14X, 53HTHICKNFSS RATIO ( D 
50UBLEF THICKNESS / SKIN THICKNESS)// 

6 1 5 X t 75H1.0 0.9 0.8 0,7 0.6 0.5 0.4 0 

7.3 0.2 0.1/) 

00 50 J = It 3 8 
00 40 L = 1# 10 
RSIGAV(L) * Sir.AV(JtL) 

40 CONTINUE 

WPITF ( 6 1 6 ) OL ( J ) t BSIGAV 

6 FORMAT ( 1H f F 8 . 2 , Fll.5 t 9F8.5) 

50 CONTINUE 

WRITE ( 6 t 7 ) 

7 FORMAT (lHl t 17HRRQGRAM CP MP L c TFO) 

STOP 

8 WRITE (6,9) 

o FORMAT ( l HI t 19H0I VERGE NT ITFRATTON) 

CTOP 

FMD 



SINCLF-tAP ADHESl VE-BUNDEO DOUBLERS 
NON-D MENS IONAL I 7ED EU8MULATJ0N EuR STIFFNESS IMBALANCE 
BFNniNG STIFFNESS PARAMETER * 0.7 


DOUBLER AVERAGE AOHEREND STRESS / MAXIMUM AOHEREND STRESS 

PARAMETER THICKNESS PATIO IDCVJRLEP THICKNESS / SKIN THICKNFSS) 



l.C 

0.9 

C .8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.0 

0.18919 

0.20588 

0.22581 

0 .2 500: 

0.28CC9 

0.31813 

C. 36342 

0.43 750 

0.10 

0. 1932? 

0,210^6 

0.23082 

0.25565 

C . 28642 

0.32554 

0. 37639 

0 *44 7 2 3 

0.15 

0 . t 9526 

C .21263 

C.233’6 

0 . 25351 

0.2 896 J 

0.32925 

0.38115 

C • 4 5 ?C 9 

0.2C 

0. 19732 

C . ’1 49 2 

0.27592 

0 .26149 

L. 29294 

0. 33299 

0. 33543 

C . 45 o95 

o,25 

0.19941 

C .2 ! 723 

0.2385C 

0.2643 1 

C. 79524 

C . 33674 

0.33971 

C • 46 161 

0.3" 

0.20151 

0 . 2 1957 

".24111 

7 « 2 6 7 2 4 

r . ?9954 

0.34051 

0.3940 1 

0 • 46 6 bo 

0. 35 

0.2036 3 

C . 22 192 

0 .24374 

v . 279 19 

C.3 f ?K 

: . 3443 1 

0.39832 

0.47150 

0.60 

0.2057 7 

C . 22430 

C .24637 

0.2731 7 

C . 30626 

C . 34312 

r • 40 ? 6 4 

0.47633 

0.50 

0.21011 

0.22911 

C. 25176 

?. 27918 

7.3 1 30 4 

0. 35 57 8 

: .<*1 129 

0.48595 

0.60 

0.21656 

3.234L t 

0 .2 572 1 

' .28 52 9 

r. 3199 . 

0.36351 

0.41997 

0.49552 

0.70 

0.219C4 

0.2 38 9 9 

C . 26275 

0.29143 

C . 32634 

0.37 k 3C 

0 . 4’8b6 

C ■ 5 0 5 0 3 

0.80 

0.22365 

C . 2 4 4 9 

■ . 2 6 e 4 1 

0.2977^ 

" . 3 3 ? 3 9 

). 37918 

0.43741 

0.51446 

0.90 

0.22831 

C .24924 

C . ? 7 4 1 3 

0 - 304 15 

0 . 34C93 

0. 397^0 

0.4461 1 

0.52330 

1.00 

0 .2 33C6 

3 . 2 544 8 

0. ’79"? 

" . 11 C 5 3 

C . 34814 

0. 395CG 

C . 4 5 4 8 C 

C . 5 3 3 1 1 

l . 50 

0.25795 

C . 2*8 1 8 ? 

: . 3 1 c : 3 

0 . 34 3 7 7 

0 . 38466 

0.43494 

0.49777 

C . 5 7 76 3 

2.00 

0.23465 

: .3109? 

0 . 34174 

i .37324 

0.42190 

r . 4 7 4 6 9 

C.5392 0 

0.61 380 

2.50 

0. 3 1 2 3 ? 

0 . 34 144 

- .37408 

0.41339 

0.46915 

0.51346 

0.57033 

0 . 65 '> 18 

3.00 

0 . 34742 

0 . 3 7 3 ■ ; i 

C . ur nC 9 

C . 44 36 1 

ij . 4 fj “> 4 -i 

0.55056- 

C . 61 46 0 

C • 6 6 9 7 C 

3.50 

0. 37""7 

C. -VC 510 

1.44165 

0 .43^29 

j . 6 3 ” v (3 

0.53552 

0.64302 

0 . 71 9b 1 

6.09 

■0 *4 rt 309 

" .43725 

:.4?4?^ 

7.516 14 

C . ‘ f>4 5 » 

0.61 309 

0.678 *C 

0 • 74590 

6.0" 

0.46535 

•3 . 49904 

5 .5 3 7^4 

0 . 57967 

0.62527 

".97548 

C .73 02 0 

0 . 7 u 9 7 4 

6.00 

0 • 5 2 3 C 9 

0.55827 

r .5053 1 

0.63492 

0.6775! 

'0 . 7?319 

" . 7 7 L 9 2 

0 .82 400 

7.00 

0 . 57775 

* . r> 1066 

0 . 6455 1 

0.63235 

0.7213? 

0.76240 

0 . R 0 5 5 6 

0.65094 

8.00 

0.6257? 

0 .66^59 

0 . 68835 

0.72263 

C . 757<j 7 

0.79451 

0.5327 1 

0 . 6 7 2 3 7 

9.00 

0 .66771 

0.69630 

. . 725 U 

" .7 5oS0 

c . r* R ? '* 

C.8? 102 

0.3546 1 

0.38953 

10. oc 

0 . 70424 

0 .73049 

C . 75749 

: . 785 17 

0 . 8 1370 

0. 84296 

C . * 729 5 

C- Q 036C 

15.00 

C . 3 2 6R 5 

0 . 34 1 5 9 

C .b5 35 .‘ 

0.87551 

C. 892 B 1 

0.91027 

0.92739 

0.94560 

20 .00 

C . 8 8 66 9 

: . 89 723 

C.9C 34 ’ 

0 .9196 9 

C . 93 100 

9<.?3C 

C . c 5 3 7 5 

C .96525 

25. OC 

0.92 042 

r , 9 28 ?° 

C."?6CS 

C • 9 4 4 C 2 

0.95196 

3.95992 

" . “6790 

0.97590 

30.00 

0.94133 

0.9^716 

C. 95^00 

0.95884 

0. 9 64 7iO 

J. 97956 

0.97648 

C * 9 8 2 2 1 

35.00 

0.955C4 

^ , 96 r 52 

r .96400 

3.96843 

7.4 72 97 

0 . 97 746 

0.96183 

0.98641 

60 • 0 " 

C .96446 

, 96 8 ~.r 

0.97155 

j. 97 r 0 9 

0 • 97564 

C. 932 12 

0.935 7l 

0.98928 

50 . 00 

0.97620 

6.97357 

" . 9 R C 9 5 

0.98333 

0.98562 

C . 0 3 8 C 2 

",99"43 

C .99283 

70.00 

0 .98720 

C . 9 i 3 4 3 

C .9896 7 

■y . 90 ,'' 0 7 

0.99227 

0. 99 356 

0.99486 

0. n 9615 

1 00.00 

C .99343 

0 ♦ 9 9 4 9 

C .99475 

0.99541 

0 . 9960 7 

C. 99672 

0.99738 

C • 9 9 8 C 4 

2 00. CO 

0.99823 

0. 783 46 

•* . 00^6? 

0.993KO 

0.99897 

0. 79914 

C .99931 

0 . 99949 

500.00 

0.9997? 

0 . 99975 

0.99977 

0 , 99 93C 

C .9993 3 

0.99986 

9 5 9 

0.9999 2 

1000 .00 

0 .9999 3 

v, • 99994 

U. 90994 

r. . 999 95 

0. 99«4b 

0 . 99996 

0.99959 

0 . 9999 2 


0.2 0.1 

C. 53846 0.70000 
0.54929 0.71030 
0.55465 0.7152b 
0.55997 0.72014 
0.56525 0.72501 
2.57048 0.72971 
0. 57563 0.73633 
C . 58 :R2 C. 73885 
0.59C97 0.76760 
C.60C93 0.75601 
0.61069 C. 76606 
0.62G25 0.77178 
C. 62959 0. 779L8 
0.63872 0.73627 
0.6811? C. 81766 
0.71811 0.86275 
0.75009 0.66337 
C. 77762 0.88032 
0.80128 0.89638 
0.82166 C. 90616 
0.85668 0.92668 
0.37933 0.93797 
C. 39860 0.96815 
0.91363 0.95602 
C. 92539 0.96222 
0.93506 0.967 L 9 
0.9636C C. 93173 
0.9767B 0.98831 
0.98382 0.99192 
C. 98315 0.99608 
0.99095 0.99568 
C. 99286 0.99663 
0.99522 0.99761 
9.99763 0.99372 
0.99369 0.99935 
0.99966 0.99983 
0.99996 0.99997 
0.99996 0.99997 


SINGLE-LAP AOHES 1 VE-BDNOEU DOUBLERS 
N3N-DI MENSIUNAL I ?E0 FORMULATION FOR STIFFNESS IMBALANCE 
BENDING STIFFNESS PARAMETER = 1.0 


DOUBLER AVERAGE AOHEREND STRESS / MAXIMUM ADHEPENO STRESS 

PARAMETER THICKNESS RATIO I DOUBLER THICKNESS / SKIN THICKNESS) 



1 .0 

C ,9 

0.8 

C. 7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

O.C 

O.25CC0 

0 ,21 '21 

C . 29 4 L2 

0.32268 

0.36714 

0.43000 

0.45455 

0.52632 

0.62500 

0.76923 

0.10 

0.25472 

C .27543 

C. 29978 

0.32882 

C. 3fc4"3 

0.40 760 

0.46790 

0 .53532 

C • 6 3 4 1 6 

0 « 7 7 6 81 

0.15 

0.257U 

: . ? 7«( 4 

0 . 30 2 64 

0.33196 

0.36769 

4ii 42 

C . 46 70 7 

w. 53978 

C .63862 

0.78051 

0.20 

0.2595L 

C . 23066 

0.30851 

0.33^1? 

0. 3 70"6 

0.41523 

C . 4 7 1 7 3 

C .54423 

0 . 64306 

0.78412 

C . 2 5 

0.26193 

^ . 283 30 

0.3C34C 

0 . 33629 

0.37445 

0.41 106 

r .4 75 30 

0.54365 

0.64753 

0.78765 

0.30 

0 . 26436 

^ 

0.3113" 

C. 34147 

n . ? 7 T94 

C . 4 2 2 3 3 

47954 

0.55304 

C . 65 1 38 

0 . 791 1C 

0.35 

0 .25632 

0.23363 

C . 31422 

C . 34467 

0. 33145 

C .426 71 

0.49363 

0.55741 

0.65618 

0 . 79449 

0.40 

C .2692 3 

?n 31 

C .31 716 

0.34738 

0. 33496 

C. 43054 

0*48782 

0.56175 

0*66 043 

C . 79779 

0.50 

0 .2/426 

C . ’ ? 6 7 3 

0 . 32307 

0. 354 3 3 

39 2-: - 

C • 4 3 6 2 J 

0.49635 

0.57034 

0 . 66376 

C . 30425 

0.60 

0.27931 

c . 3 : 2 7 1 

0 . 3 29., 3 

:. 36CO J 

0 . 390:3 

C . 4 4 5 9 6 

C. 5342^ 

0.57802 

0.67638 

C . 3 l 039 

0.70 

Q .23441 

‘.33775 

C. 335C6 

7 . 36788 

0. 4Cnl 8 

3.45352 

0.51237 

0.58717 

0.6843C 

G • 9 1 6 27 

0.80 

0.23953 

0.31375 

j , ? 411’ 

0.37396 

0.41330 

C .461 17 

0.52044 

0*59540 

0.69231 

0.82191 

0.9C 

0 . 2 943 : 

0 , 3 1 Q : C 

".34724 

0.38058 

0. 42C 44 

C • 46579 

C. 62945 

0 .6C3 49 

0 . 70000 

0*92731 

1 .CO 

0.30" l 1 

0 • 3>474 

, . 35344 

0.3B727 

C . 42 758 

0.47639 

C. 5 1633 

C .6 1 145 

0.70710 

C. 83248 

1.50 

0.3?73’ 

0.35397 

: . 3 34 7 8 

3.42031 

0.46328 

0.51337 

0.57432 

0.64912 

C. 74031 

C . 855 3 l 

2.00 

0.35560 

'.33403 

C .41676 

0.45453 

C .49843 

0,54994 

0.61077 

0.60312 

0.76972 

0 . 97395 

2.50 

0. 33465 

C.41467 

0 • 4 4 h H C 

0 .‘♦0 780 

0. 53257 

3.58420 

0.64396 

0.71351 

0.79460 

0 . 99929 

3.00 

9.4141 1 

0.44534 

0.48049 

: .520 20 

0.56516 

0.61622 

0.67428 

C . 7495C 

0 . 0 1 60 C 

0. 90 2 04 

3.50 

0.44353 

0.47568 

".51143 

0.55134 

C • 5 9 55 6 

3. 64 8 03 

3. 70177 

0.76439 

0.33438 

0.91272 

4. CO 

0.47275 

: . 5: 534 

7 .54129 

0 . SB" 96 

0.62478 

0 .07 312 

3.72657 

0 .78551 

0.85035 

C . 9 2 1 7 6 

5.0C 

0 • 5 ’ 89 7 

0 . 56! 86 

C , 59679 

0 . 63497 

0.67624 

0. 72C34 

0 . 76 3 3 7 

0. 82068 

C. 37637 

0 . 93607 

6.0" 

0.56036 

: . 5 1 2 3 6 

C . 64594 

" .68178 

0 . 71963 

0. 76 016 

u. Q 0 30 1 

0.34836 

0.09621 

C. 94679 

7 0 

0.62743 

0.o5730 

".68872 

0.7217? 

0 . 75636 

C. 79261 

0 . 8 1C 6 2 

C.87C15 

0.91177 

C. 95502 

8.00 

3.66860 

0 .49644 

0. 72542 

0 . 75557 

0. 73684 

0.91937 

C. 85310 

0 .38 794 

C. 924 L2 

C . 961 47 

9.00 

0,70455 

: . 73-25 

0.75680 

0 . 7 H 4 1 4 

0.31243 

C. 94156 

0.*71 55 

C. 90234 

0. 93406 

0 ■ 9666 l 

10.00 

0.73577 

c .75937 

0.73354 

C.8CC43 

0.33394 

0.66CC7 

0 .80675 

0.91417 

0.94218 

0.97071 

15. OC 

0.83997 

0.35625 

7 .8706? 

0.38625 

G.9C204 

0. Q l 799 

0. 93409 

C. 95028 

C.9667? 

0.99329 

20.0C 

0 . 89419 

0.97453 

" .0 1492 

0 .925? 7 

C. 93579 

C .94c»33 

0. 95^C? 

0.96769 

0.97841 

0.93913 

25.00 

0.92521 

0.93251 

0. 93994 

0. 94738 

0.95484 

0.96232 

0.96932 

0.97733 

0.98479 

0.99240 

30. CO 

0 .94440 

0. J 4 9 9 2 

7.95544 

0 . 96C98 

0.96653 

3.972C8 

0.97765 

C.9R314 

0.98877 

0.99439 

35. CO 

3 . 0 5 7 1 1 

C. 961 38 

0.96565 

: . 96993 

0.97421 

C. °7349 

C. 99271 

C.987C4 

0.9? 136 

0. 99568 

40 .00 

0.96592 

C. 9t,9?2 

0.97272 

j. 97612 

C. 97943 

0.93236 

O.Q0629 

C .98972 

0*99315 

0.99658 

5C.00 

0. 97703 

0.17930 

" • 98 l 6t 

0.98 389 

C.986 1 1 

C. 93843 

0.99C 75 

0.99307 

0*99533 

0.99769 

70.00 

0.93752 

0. 93667 

J .949 94 

C . 99 120 

C. 59247 

0.99373 

0.99499 

C .99624 

0 . 99750 

0. 99875 

100.0C 

0.99355 

0 .094 ?" 

3. 994H5 

0.59549 

C. 99b 14 

0.99673 

0.99743 

C . 99 0"* 7 

C. 99 372 

0.99936 

200.0C 

0 • 99 3 3C 

0.99847 

0.99064 

C . 9938 1 

C. 99898 

C. 99915 

0 .99932 

0.99949 

0.99966 

0.99983 

500.00 

0 .9997? 

7 .99975 

0.99977 

"■ . 99930 

C. 9 99 13 

C. 99986 

C . 99939 

C .99992 

0.99994 

0 . 09997 

10CC.00 

0.99993 

C . 09994 

- . 99994 

3.99^95 

0.99996 

C. 99996 

0.99989 

C .99992 

0 . 99994 

0 . 99997 



S1NGLF-LAP ADHF SIV6-«JNDfeO DOUBLE R S 
NON-DIMENSIONAL 1760 FORMULA! 1 FOP STIFFNESS IMBALANCE 
RENDING STIFFNESS PAnAMFTEK * 2.0 


DOUBLER 

L/T 

PARAMETER 


4VFPAGF AOHFRENO STRESS / MAXIMUM ADHERENO STRESS 
THICKNESS RATIO (OGURLFR THICKNESS / SKIN THICKNESS) 


0.0 
0.10 
0 • X 5 

0.20 
0.2 5 
C • 30 
0.35 
0 • AO 
0.50 
0.6C 
0.70 
0.80 
0.00 
1.00 
1.5C 

2. CO 

2.50 

3. CO 

3.50 
A. 00 

5.00 
6 . C 0 

7.00 

8.00 

9.00 

10. OC 
15. OC 
20. OC 

25.00 

30.00 

35.00 
AO. CO 

50.00 

70.00 
IOC. 00 
200.00 
500.00 

1000.00 


i.C 

• AOC 0 
, AC 53 A 
.ACM jA 

. AIC73 
.413^6 

• A I 5 1 6 
. A 1836 
. A 2 l 5 7 
.42693 
. 4324} 
.43732 

• A A 3 2 A 
, A A 6 5 A 
. A5AC 7 
•A810? 
. 5C 7 5 4 
. 5 3341 

• 5 5 B*»A 
.58249 
. 6354 5 

• 6 A 79 1 
.63556 
.71859 

• 7 A 7 5 1 
.77273 

.79470 
. 85 q 7 B 

.9 1 r-i i 

.93537 
.951*8 
.96170 
.96921 
.97*06 
.95318 
. 99334 
.99034 
.99972 
.99993 


0.9 

. 42553 
.4311 A 
. 43397 
.4367? 

. A 396*. 

.44? 47 
.44529 
.44811 
. 45376 
.45939 
. A65 r 2 
. A 7 C 6 A 
.4 76 25 
.43184 
. 5 0951 

. 53 54R 

.56254 
.5975? 
.41131 
. 4 3 * F 3 
. 6 7 5 c ; 
. 7 1 1 5 P 
. 7 A 2 2 
.7 694 7 
. 79297 
.51333 
.83224 
.91961 
.94 174 

.95593 

.96551 
. 97 228 
. 9^097 
.98937 
. 9 >4^6 
. 94 8 50 
.9^975 
. 9 >904 


c.a 

Q. 45455 
A.^oOA 1 
0.46333 
0.46632 
C .46925 
v. 4 7224 
0.47518 
C. 478 12 
0. 4*398 
i.4093? 
0.49564 
A . 591 44 
0 .5C 721 
0.51298 
C .541 19 
" , 5 (jB 4 C 

: . 59442 

0 . 6 1 9 1 C 
' . 6423 7 
C .66421 
C. 7C36C 
C. 73767 
C . 76697 

n . 792 1 3 
C .81375 
r . 8323* 
** . 0968 2 
0.92326 
0.94814 
0.96079 
C • 9 6 0 3 2 
0.97535 
C . 0 8 3 08 
•; .99 056 
C . °95C3 
C. 99867 
0*99978 
0.99994 


C.7 

0.43781 C 
C. 49393 C 
0.49701 C 
0.50CC6 0 
0.5031 1 r > 
0.5C621 0 
J. 509 ? 8 C 
0.51229 0 
C. 51 83 A 0 
0.52435 0 
C. 5303 3 0 
0.53627 j 
C. 54216 0 
0 .54801 : 
:. 57656 0 
C .60 37 2 C 
0.62937 C 
0.65344 C 
0.6758V C 
0.69675 0 
0.73391 0 
0.76561 0 
0.79256 0 
C . 3155C 0 
0.83507 0 
0.R51S5 C 
0.90753 0 
C. 93 70 £ : 
C. 95455 C 
C. 96566 C 
0.97314 0 
0.97832 C 
0.98 5 1C C 
0.991 7 A C 
0.99569 l 
0 .99884 C 
99981 


0 . 99095 


0.6 

.52632 0. 
.53265 C. 
.53583 ‘ 
.53878 
.54212 
.54525 
. 54843 
.55154 
. 55772 
.56385 
.56992 
.57593 
. 581P9 0 
. 5877 7 C 
.61623 0 
.64293 0 
.66779 0 
.69 08 3 0 
.71 2C 7 0 
.73161 0 
. 76597 0 
. 79485 c 
. ft 1 9 1 2 0 
. 83969 A 
.85 695 C 
.87166 C 
.92^3 6 C 
,64694 C 

.96099 C 
.97053 0 
97696 0 
°P 14? C 
,9tt724 0 

90293 C 
996 31 C 
99900 C 
99983 0 

99996 0 


C. 5 

57143 0, 
5779: 0, 

5 H 11 4 0, 
.58435 C 
,58753 0, 
,59070 C 
,59385 0 
. 59 7C 5 C 
, 6 r 3 26 C 
,60940 l 
.6 1 546 0 
.62144 : 

.6 >734 0 
.633 15 0 
.66096 0 
,6B6'>0 0 
.71011 C 
.73160 j 
.75117 0 
. 76 ft 9 7 C 
. 7 99 86 0 
.8254? 0 
, 84 66 7 0 
. 96435 0 
. 8793 1 0 
.8920 0 C 
. ->3324 0 
.954R4 C 
.Q6744 0 
.97542 C 
. 9337u 0 
.9845? C 
.98937 0 
.9 44 11 * 
. 99693 C 
.99917 C 
.99936 C 
.99996 C 


C • 4 

. 62 500 
.6314 7 
.63469 
. 63 73 7 
.64 102 
. 64415 
. 6 4 72 6 
.6 50 34 
.65648 
. 65247 
.66335 
. 67412 
.67979 
.433 <6 

.71157 
.7354 1 
. 75683 
. 7761 1 
.79336 
.8^395 
, 53567 
, 85730 
,3751 6 
,38996 
,9r?33 

,91276 
, 94636 
, 96379 
,97391 
,98032 
,98456 
,0876? 
, 49 150 
.99529 
,97754 
,99934 
. 99939 
, 99939 


0.3 

. 68966 
.69584 
.69891 
.70193 
.70491 
.70 736 

• 7 l C 7 o 
.71367 
. 7 L 93 9 
.72475 
.73037 
. 73566 
. 74 084 
. 74588 
.76934 
. 79010 
. 6 C 3 4 7 
.32471 
.8391 1 
.35190 
. 67343 
. 69 074 
.90479 
.91632 
.92590 
.93393 
.95960 
.97277 
.98041 
.96517 
. 98842 

• 9907 1 
.99 363 
.99447 
. 99« 16 
. 99950 
.99992 
.99072 


0.2 


0.76923 0. 
0.77463 0. 
0.77729 0. 
0.77990 0. 
0.78246 0. 
0 . 7 8499 C. 
0.73748 0. 
0.73993 C. 
0.79469 0. 
0.79934 0. 
0.80393 0. 
3. 398 31 0. 
0.81256 0. 
0.31668 C. 
C. 33550 0. 
0 . *5173 0. 
0.86576 0. 
0.37796 0 . 
0.33860 ‘ 
3.39793 
0.91341 
0.92564 
0.93546 
C .94345 
0. 95003 
0.95553 
C .97295 
0.93160 
0,985*7 
n • 990 1 1 
C. 99228 
C. 99381 
C .99576 
0.99765 
0.99877 
C ,99967 
C. 99994 
0.99994 


0.1 

86957 
37 322 
87500 
37673 
37843 
33010 
33173 
38334 
33643 
33941 
89228 
89504 
09769 
90025 
91172 
92128 
92935 
9 3622 
942 12 
94722 
95556 
96204 
96710 
97127 
97468 
97752 
93633 
99C85 
99344 
99506 
99614 
99691 
99733 
99883 
, 99939 
,99984 
, 99997 
, 99997 
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r SINGLE-LAP ADHFC IVE-RONOED DOUR L. FRS 

C ADHEREND PLASTIC HINGE F OR V A T I ON (HR foactijpF) - UNRALAwrfD STIFFNESS 
C ADHFPEfc" PENDING A*'0 c X T C N S T ON AL STIFFNESSES UNCOUPLED FOP COMPOSITES 
C NON-DIMERS IONAL 17 FP FORMULATION 

DIMENSION 01(38), S I G A V < 3 8 , 1 0 ) , FSIGAVUOI, BENDKC l 3 ) 
f S FT OVERLAP ARRAY 
OKU = 0. 

OK 2 ) = 0.1 

OK 3) = 0.15 

01(41 = 0.? 

01(5) = 0.25 
OK 6) = 0.3 
0L( 7 ) = 0.35 

01(8) = 0.4 

PL ( 9 ) = 0.5 

OK 10) = 0.6 
OL ( 11) = 0.7 
OKI?) = 0.8 
OK 13) = 0.9 
OK 14) = 1 . 

OK 15) = 1.5 
n L( 16) = 2. 

0L( 17) = 2.5 
OK 18 1 = 3. 

OK 19 ) = 3.5 
OK 20 ) - 4. 

0L( 21 ) = 5. 

OL ( 22) = 6. 

OL ( 23) = 7. 

OL(24) = 8. 

OL ( 25 ) = 9. 

0 L f 2 6) = 10. 

OL ( 2 7 ) = 15. 

n L ( 2 8) = 20. 

OL ( 29 ) ^ 25. 

01(30) = 30. 

0L( 31) = 35. 

OL ( 3 2 ) = 40. 

0L( 33) = 50. 

01(34) = 70. 

0L( 35 ) = 100. 

OL(36) = 200. 

01(37) = 500. 

01(38) = 1000. 

C SET BENDING SMFFNFSS PA D A M c T r P ARP AY 
RFMHK ( L ) = 0.1 
8 E N n k ( 2 ) = 0.15 
RENnK(3> = 0.2 
8ENDK(4) = 0.3 
8 FNP K ( 5 ) = 0.5 
8 E N n K ( 6 ) = 0.7 

8 E N D K ( 7 ) = 1. 

8 f n n K ( 8 ) = 1.8 

qFMnK(o) = p , 

8ENOK ( 10 > = 3. 

8FN0K Cll) = 5. 

«ENOK (12) = 7. 

8Fmp,k(13) = 10. 

C SELECT PENDING STIF«=N C SS P^/mftpd 
00 50 m = l, 13 

C SET A DH ^ c ENO THI r KH c SS RATIO 
TP - ll. 

90 30 K = 1 , 10 

TR = TP - l. 

T in VT 4 = TO / 10 . 

D10V94 = T 1 0 v/T 4 . ** 3 

C SET I K T T A t ESTIMATE AV c ° AGF AOHFP c Nn STRESS 
V 3 = l . / (1. 4 (3. / 8CN|pK(M>) * T 1 n V T 4 ) 

S I G 4 V < 1 , K ) = V 3 

90 20 J = 2, 38 

r SET NON-OI MENS I ON A L f 7 C D OVERLAP 
V2 = °L ( J ) ** 2 
N = 0 

L = N + 1 

IF ( N . ot # 500) G9 T9 a 

c pPOGPAM i)SF^ L(npH8l8D> / T(SKTN) AF l/t PA T IO 
V1Q = V3 * V? / B FN n K ( *) 

C TO USE L/t patio FOR no»)3L E° TMS T F4n pF SKIN, US r 
C V10 = V3 * V 2 * D l n V 04 / REMPK(v) 

r ppQ c i vplY-SHP POP r E9 t^NOS OF SK IN- POMP L F 5 C nM 8lNATION 
V8 = l. SORT { V 1 0 ) + V10 / 6. 

FOR BUILT-IN p NOS HF SK I N- DO! I* L F* COMBINATION 
V8 = 1 . + (SOPT(viO)) / 2. + V 1 9 / 24. 

EVALUATE ECC p NTP K tty OAPAtrrpp 
FK = 1. / VR 

PF-ESTIM4TF A V c R A G r Stbf^p 

V4 = l. / (1. ♦ (3. / 8 r V 9 K ( m ) ) * (T10VT4 / V3)) 

r CHEEK CONVERGENCE OF AV r PAnr STRESS 
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p = V3 / V4 

Jp ( (1. 0001 .C T . R) .AND. (0.9990 ,LT. 91 ) GO TO 20 

V 3 = V4 
GO Tn 1 

EVALUATE AVFRAGF SK!»I STRFSS AT LIMIT LOAD 

( V T FLO STPFMGT H po° MFT4L< . ULTIMATE STRENGTH FOR COMPOSITES) 

20 STOAV (J,K) = V3 
30 CONTINUE 
C PPJMT hut T 4BUI ATT<->v|c 

WRITE (6,5) BFNOK(M) 

5 FORMAT ( 1H 1 , 51/ ), 30 X , 35HS I NGL E-L A° AOH c S I VF-BONOEO 0ni)BL c P5// 

1 20X , 55HVnN- , iT'' c NSinN4LI2Pn F JR MU L A T I ON FOR STIFFNFSS ImoalANCF// 

2 30 X , 30HBFN01NG STIFP\’ESS PARAMftER = , F5.1 /// 

3 3X , 7H0niJPLFR. 18X, 49HAVEP.'0P A OH E RENO STPPSS / MAXIMUM AOHERPNin 
A STRESS/, 5X , 3HL/T / , 3X, OhpaRAMETPP, 14X , 53HTHI C KNF SS RAtjo (D 
50UBI PR THICKNESS / SKIN THICKNESS)// 

A I5X , 75HL.0 0.9 O.B 0.7 0.6 0.5 0.4 0 

7.3 0.2 0.1/) 

no SO J = 1, 33 

DO 40 L = L, 10 
PSIGAV(L) = S I C A v < J , L I 
40 C n NT I M( I F 

WRITE (6,6) OUJ), R S I G A V 

6 rnpriAT ( LH , F3.2, FU.5, 9P8.5I 
50 CONTINUE 

M'»I TC (6,7) 

7 FORMAT ( LH 1 , 17HPPDGP4M COMPl'Tcn) 

STOP 

3 W»ITc (6,9) 

9 FORMAT (IMi, 19H01VPRGENT IT c RATIOM| 

STOP 

FNO 



SINM.6-LAP AUHES I Vt-0UNDtD DOUBLES 

VIN-DI 'ENSIDNAL I 7F9 FORMULATION FOR SUFFNFSS IMBALANCE 
TENDING STIFFNESS P AKAMC TER - 0.7 

AVERAGE AOHERfNO STRESS / MAXIMUM ADHr*ENO STRESS 
THICKNESS RATIO (DOUBLER THICKNESS / SKIN THICK NFSS) 


0. 1891* 
0 . 1 * 7 ? 2 
3.2C15 1 
3.2C577 
0.21011 
0.21454 
0.21*34 
0.223^5 
0.23 BOS 
0.24273 
0 . 25261 
0.26313 
0.273 77 
0.20466 
0 . 36251 
0 .60366 
0.46535 
0. 5230«> 
0.57775 
0.6257? 
0 .7C627 
0.763 r 6 
3 .807M* 
0. 86C37 
0 .86606 
0 .8351 1 
0.041 37 
C .96646 
0.97620 
0.9 8 2**6 

0.98716 
9.09001 
9.99 346 
0.99657 
*.9*8?0 
0. 999 6 6 
9 .9999 3 
0.99993 


I?1 49? 6, 
. > 1 n 6 7 ° 
‘,>24^0 0 
. 2 2 9 ] 1 0 , 
'.>3401 r 

t ? pm .* 

' . ? S, 4 0 > 0 
1.75447 •: 
^.>6513 : 
.77917 : 

. > 9 7 4 J 0 
’ .2 *9.; 9 0 

.11,9? 0 
). 1 73: ' 3 

; . 4 3 7 > 1 : 

.499 54 . 

> , 5 } Cl ? 7 ' 

; .->l"6 5 «' 

: . 6 5 6 ci : 
■’.73:44 ' 

‘ . 7 3 4 r ’/ r 
. H ? 6 2 r 0 

: . 66*73 c 

■ ,37^C5 C 
.89725 C 
: . >4 7 19 c 
"',94709 c 
j. 07E 57 C 
-.^■1457 r 
0 • 9 8 o 4 4 0 

3.99i:i - 

:.9Q4io 0 
C . 994* 1 0 
3.99346 : 
' . 9996 L F 
G. 99 994 C 
:.*99*4 : 


,22581 
, 2 3 5 ^ l 
,24111 
,24639 
,26176 
.25 721 
,26275 

,269/,l 

,279*2 

.29174 

,3:336 
. 31624 
.323^0 
.34175 
. 4CH: r 
. 47476 
. 507 c 4 
.*9*^21 
.64661 
. 6 0 8 2 r 
. 7 4 766 
.8072-3 
. 84 3 8 l 
.H 7122 
. 8 ‘ * 2 1 5 
.9336 
. 95295 
. 9715 h 
. 980 9 5 
. 98 o? 9 
,98972 
.9921 1 
. Q9 476 
* 9 Q 7? 6 
. 09 *? 6 5 

.99994 
.9 *9 94 


: . ?s:r : o, 

7.26139 G 
".26 72 3 0 
3.27316 0 
0.27913 C 
0 . 2 4 5 2 9 : 
:.?9i4b ■" 
0.2 977 9 ■; 
".31058 3 
'.3>2h6 0 
3.33703 0 
'. 3606 4 C 
0. a o4>4 C 
_ . 37 00 5 0 
0.44 8 f j4 6 
: . 5 1 6 9 6 3 
:.5796 7 0 
0.63492 r 

,6B?8 3 c- 

".7226’ C 
3.70522 c 
~ . 8 7 q 9 7 0 
C • 362 r 7 0 
*' . 886 86 0 
0.905*4 " 
A . 919*4 C 
3.94 8H1 r - 
0.97508 0 

G • c H 3 3 3 0 
C.8H0C1 C 
0. 9*1 0 1 C 
■:• 99301 0 
0 . 9 9 5 4 ? 
r, 9976C 0 
, # qo 8 3C 0 
0.99949 0 
"-.9898 5 C 

0.99995 f 


2*:: c 

,29294 0 
,29855 0 

,32 62 5 ■* 

, 31 3C4 : 

, 318): c 
, 32404 0 
,33339 0 
, 34 *1 3 2 
,36262 0 

.37725 0 

,39?05 r 

,40694 : 

. 4 2 i 7 2 C 
, 4°5 7 L 0 
,54453 ; 
.62537 0 
.67 751 j 
. 72 1 32 0 
. 75 7b7 0 
.81 373 C 
.06308 C 
. 0HIS5 r 
. 8 r 2 fc 3 0 
.91055 0 
.93° *7 : 
.90465 3 
.078*4 0 
.9856? C 
.00473 C 
. 9° 2 3 0 0 
.9940 1 0 
. 0 9 o C 7 0 

.09 794 C 

.99897 C 

# 999T4. : 

.09906 : 

.99900 : 


,31818 : 
.3 329 3 ° 
, 34 C 5 l C 
.34811 0 
.56570 0 
.30351 0 

,3713c 0 
.37910 C 

.3980} 0 
.4 1 Co 3 C 
. 42 6 92 C 
.44291 ; 
.45884 0 
.47460 C 
.55:5? C 
.41003 C 
.67540 0 
. 72 3 19 0 
. 7o?4G 0 
.79 4 51 0 
.0428? C 
.6/862 C 
.9: 375 0 
.91 044 C 
.33194 T 
. 942 J 5 C 
.97055 C 
. 9H ? 2C C 
. 93 8C 2 C 
.99145 : 
.97350 ( 
.975 Cl C 
.99673 C 
.99029 : 
.99014 ; 

.99970 : 

.99996 I 

.09996 ^ 


.34042 , 
. 30 54? i 
.374:1 1 
. 4026 3 1 
.4 1129 
.41997 
. 42 06 4 
. 43741 
.45400 
. 472 1 1 
.43*26 
.50619 
. 4?>0 5 
. 53910 
. 61456 
.67329 

.73:2? 
.77192 
.80556 
.03271 
.07293 
.90085 
.92010 
. 9 346 3 
. 94540 
,9537b 
.97642 
.90575 
.99343 
.99316 . 
.9940 7 
. 9940 l 
j . 9 9 7 3 0 
' .9906 3 
.98931 
' .98902 
; . 9 799 7 
; , 88 99 7 


. 43 750 C 
.45694 0 
. 46 ofe5 0 
.47632 C 
.48595 0 
.49552 C 
.50 50 3 0 
.51444 0 
.53310 9 

.55129 0 

. 56 70 0 0 
.58619 0 
.60 20 1 C 
. 6 L 0 0 3 C 
.6 857 4 0 
.7459 3 0 
.78974 C 
,0 240-C C 

.05094 0 

.87237 0 
. 9C 3 6 5 0 
.924 70 0 
. 93 *7 8 C 
.9507? 0 
.95093 C 
.96526 C 
.98231 0 
.98924 C 
.98233 0 
.89487 0 
.99615 0 
(.99701 0 
J. 99 80 4 0 
99R97 C 
.99^49 0 

‘.99037 0 
.99998 C 
1 . 99 99 R C 


,53846 0. 
.55996 G. 
. 57C 40 0. 
.50082 0. 
. 59C 97 C. 
.60093 0. 
.61069 0. 
.62025 0. 
.63076 C- 
.45637 0* 
.67310 0. 
.60896 0* 
.7D39o 0 . 
. 7 10 1 3 0. 
. 77 759 0. 
.82165 0. 
.05440 0. 
.37933 C. 
.09840 0. 
.91343 C. 
.93503 0- 
.94951 G, 
.95965 0, 
.96 7C2 0, 
.97254 0, 
.97671 C, 
.90820 C. 
.99204 0, 
.995 22 0, 
.99653 0, 
.99744 0, 
.99901 0 
.99069 0, 
.99931 0, 
.99966 0, 
.99991 C, 
.99991 0 
.99991 0 


, 7C00G 
. 720 ?0 
, 7297 1 
, 73884 
,74760 
, 75601 
, 76406 
.77173 
.78625 
. 79955 
.81176 
.82293 
. 83 330 
. 94279 
. 00036 
. 906 L 1 
. 92449 
. 93797 
.94615 
.95602 
.96719 
.97458 
.97965 
.98340 
.90620 
.98835 
. 99406 
.99643 
.99761 
.99829 
.Q9872 
. 99900 
.999 35 
. 99966 
. 99983 
. 99996 
. 99996 
, 99996 


SINGLE* LAP AUHfc S I VE- w ONOFO DOUBLERS 
NON-DIMENSIONAL! 7 E H FORMULATION FDR STIFFNESS IMBALANCE 
BENDING STIFFNFSS PARAMETER = 1-0 


DOUBLER 

L/T 

PAR AMETFR 


AVERAGE ADHf REND STRESS / MAXIMUM A DHE 7 END STRESS 
THICKNESS RATIO (DOUBLER THICKNESS / SKIN THICKNESS! 


0.0 
0.10 
0.15 
0.20 
0.25 
0.30 
0.35 
0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
i .00 
1.5C 
2. CO 

2.50 

3.00 

3.50 

4.00 
5 .00 
6. SC 

7.00 

8.00 

9.00 

10.00 
15. CO 
20.00 

25.00 

30.00 

35.00 

40.00 
50 .00 
7C.OO 

100.00 
2 00 . C 0 
500.00 
1000.00 


0.25000 
0.26951 
0.26436 
0.26020 
0.27426 
0.27931 
0.23441 
0.20950 
0.3" Oil 
0. 31CS4 
0.32177 
9.33208 
0.344 16 
S.3555B 
0.4 141 > 
0.47278 
0.52097 
0.50036 
0.6? 743 
0.66860 
0.73574 
0.78628 
0 .82437 
0.85360 
0.87630 
0.09421 
0 .94443 
0.96591 
0.97 jzr 
0.90335 
9.93747 
0.99023 
C. 99357 
0.99661 
0.99330 

3.99956 

0.99993 

3.99993 


S.27C27 
0.23069 
0.20595 
' . 2 7 1 3 l 
0.2967? 
0.33221 
0 . 3 ’' 77 5 
0.313*5 
0 . ? ? 4 7 3 
0.33429 
0. 34003 
C , 3599? 
0.37104 

: . 334:6 
' . 44 5 5 5 
7 . 50 5 36 
r , 56 1 56 
0.61235 
65 73^ 
0.6*644 
0. 75935 
^ . 006'* 7 
0. 34 102 
C . 06 76 7 
: . 00 3 30 
0.93446 
3.94986 
■3. 9693 1 
0 . * 7 9 3 :• 
0.90502 
C .7 0 8 7 3 
G . Q q l 2 l 
0 *9^421 
0. 99 6-9 5 
0. 99847 
0 .99961 
,99994 
C . 99994 


0.2*412 > 
'.30551 : 

0.3113-° 9 
p . 3 1 7 1 6 :> 

I .32307 0 
32 9C 7 : 

0.3*505 ) 

3.34113 0 
0 • * 5 3 4 4 C 
C . 06509 0 
C. 3 7 047 0, 
C . 39117 3 

‘.40393 C 

'.41674 S 
C • 4 3 r 4 7 j 
0.54125 : 
0.59679 0 
C .64 59 9 0 
0.60872 ^ 
0.72 54? 0 
0.78359 C 
0.32616 0 
0 , 8573 7 0 
0.80107 0 
o . 9 ; : * 9 c 
O.9L407 ^ 

C. 95541 0 
C. 97?71 3 

0.95lt>v 0 

C • * 8 66 9 3 
r .9 099 8 0 

C . 99 > l 9 3 

0.994.0* C 
C. 99729 C 
S. 99864 0 
0.99965 0 

^.99994 ° 
0.99994 0 


3225 4 
, 3 3 61 1 
,34147 
, 34 747 
, 3^43 3 
, 36 C * 3 
, 36 7jft 
. 37 V*6 
, 30 726 
, 4CCf 4 
, 4l4t' 9 
,4>7S0 
.44106 
.45451 
.52010 
.58094 
.67497 
.68178 
.7217? 

. 7555 7 
. 40 846 
• 64 6 b 9 
. 07493 
. 89620 
.91249 
.92534 
. 960*6 
.97611 
.90389 
.90036 
.99124 
.99317 
.99550 
.99763 
.9900 L 
. 99969 
.99995 
. 99995 


0.35714 
0. 37 0*6 
C . 3 7 7 r ♦ 
0. 3b4*5 
0.390 ' J 
: . 399.. 6 
0. <*061 0 
0.41330 
0.42 7* 2 
0.441*? 
0. 4 66 1*3 
0.^703* 
0 . 4°45° 
G . 4 9 k 4 6 
0.56515 

0.62477 

0.67625 
0.71*07 
0 . 75o36 
C * 7 h 6 8 4 
0.83390 
C. 86 75 7 
G. 8921 ? 
0.91050 
0.92477 
3. 935P5 
0.96651 
0.97952 
0.98611 
0. 99CC3 
0. 99 249 
0.99414 
C.996L4 
C . 9 Q 7 9 7 
C . 99 0 * 8 
C. 99974 
0.99 ^ 9 6 
0.999*6 


0.40 CSC 
C. 41523 
C . 4 ? ? 8 8 
0 .43054 
".430>O 

0.44606 
C. 45352 
0.46117 
^.47643 
0.49154 
0.5 0 649 
: .5? 1 24 

0.53574 
0.54997 
C.ol 621 

0.67017 

0.72004 
G. 76C16 
0.7926 1 
0.81937 
0 • 8 6 0 0 5 
0.03079 
0.9395b 
0.92519 
0. 93712 
0.94542 
3.972:8 
0, 932*2 
-j. 9 3 04 3 
0,99169 
3. *9375 
0 . 9*5 12 
C . 996 7 9 
C. 90031 
0.99915 
C. 99973 
C . 99996 
0. 99996 


45455 1 
0.47123 
: . 4 7 0 5 4 
: . 4b 731 
: .49605 
: . 50 4>4 

: . c 1 2 i 7 

0 . 5 > 0 4 4 
: . 53647 
0 . 5 5 2 0 5 
0.567 34 
0.502 24 
0.5*673 
C .61030 
0.67431 
C . 72552 
'■'. 76007 
0 . 0 0 30 l 

0.33062 

C . 8 5 J 1 C 
0.00602 
0 . 9 l ") 2 7 
0. *27?9 
0 . *0992 
r . 949 5 4 
0 .*5703 
0.97764 
0 . *«6* 3 
0.99775 
0.99 036 
D. 99500 
0.9*610 
0.9974? 
0.99865 
0 .99932 
0.09983 
0 . 999*7 
0.99997 


.52632 0. 

.5447? 0. 
*55304 0 . 
.56175 0. 
.57034 0 . 
.57032 0. 
.56717 C. 
.5*54'- 0 . 

.61142 0. 
. 62697 0 . 
1.6410* 0. 
,.65623 0. 
(.66 7*7 0 . 
‘.60313 0 . 

74047 c. 
„ . 7 q 549 c. 
‘.02060 0, 
‘.84036 C. 

07035 0. 
'.8 8 794 C. 

91412 0. 

1 . 93222 C. 
;.945 19 9, 
>.95477 J< 
:.96204 r, 
:. 06770 o, 
98322 0, 
98969 C, 
?. 99307 9, 
1.995C? C 
: .09625 C , 
‘‘.99 707 C 
? . 9 9 80 7 0 
S .99399 C 
S. 99949 C 
C. 99987 0 
0.99998 0 
0.99*98 C 


,62500 
,6431? 
,65183 
,66042 
,66876 
,67683 
, 6043C 
. o9251 
.70723 
.72126 
. 7 344b 
. 7 4 690 
.75963 
.76959 
.8150* 
.95041 
.07637 
.00621 
.91177 
.*2412 
.94217 
. 95450 
.96327 
.96*73 
.97453 
.97835 
.98871 
. *9314 
.99538 
.99668 
.99750 
. *9 305 
.99672 
.9993? 
.99966 
. 99991 
. 99991 
.99991 


0.76923 
0.78415 
0.79114 
0.79784 
0.80425 
0.81039 
0.81627 
C. 32191 
0. 8324ft 
0 . 84219 
0.85113 
0 . 85936 
C . 34695 
C. «73*7 
C . *0 20 1 
C. 92174 
C. 93607 
C .94679 
C. 95502 
C. 96147 
0.97079 
0- 97709 
G. 98148 
0.98477 
0.98725 
C. 90917 
C. 99437 
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S I NC.l F-L A? AOHES [VF.-BHUDEQ DOUBLERS 
NON - D 1 MFNSIUNAL I ZED F ORMUL AT ! ON f-Cl« STIFFNESS l MRALANC6 
RENDING STIFFNESS PARAMETER = 2.0 


doubler 

L/T 

PARAMETER 


AVERAGE AOHERENO STRESS / MAXIMUM A OME RENO STRESS 
THICKNESS ° A T 1 0 JOOURLfcR THICKNESS / SKIN THICKNESS! 


0.0 

0.10 

0.15 

0.20 

0.25 

0.3C 

0.35 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

1.00 

1.50 

2.00 

2.50 

3.00 

3.50 
<>.00 

5.00 
6.0C 

7.00 

3.00 

9.00 

10.00 

15.00 
2C .00 

25.00 

30.00 

35.00 

40.00 

50.00 
7C.OO 

100. CO 
200 .00 
500.00 
1000 .00 


1.0 

AOjCC 
41?7;> 
415 16 
42 157 
4240 3 
4 3 2 42 
4 3 7 3 2 
44324 
4 5404 
4 64 8 6 
47567 
43633 
,49 72 1 
,50755 

,55847 
, 60 54 9 
,64 7<U 
, 66644 
,71359 

,74761 
,79464 
,83044 
,863 36 
,87 94 > 

.39 r : 4 

91C*7 
,95 U 4 
. 969 ? \ 
,*77386 
.9345? 
,9^?^ 
. 9 'r 7 4 
99 *36 
.9^.4 7? 
,qq? 34 
.99957 
.99991 
. 99 4 9 3 


0.9 

.42653 
.4 3^3? 
.44 >4 7 
. .4 6 11 
. 483 76 
. 4 5 ° 3 4 
. 4 4 6 r 2 
.47064 

.43 132 
. 4 9 ? 9 5 

• 5 0 4 0 4 

• 51 49 8 

.52882 
.53649 
.53755 
-.6 3'’ 

. 675C j 
. 7 1 C 9 3 
. 74222 
. 7r>447 

. d L 3 3 \ 
. 8 4 44 u 
. fl 7 J ? .0 

.3*182 
.4" 7 r 5 
, 9 1 *7 5 7 
. 95590 
.97727 
i. 9^ ;m7 
' .43 ‘. 0:7 
!. 94 4 *2 
>.99167 
; . 99 44 7 

> . o 9 7*5 
) . 4 9 5 l 
; . 9 4 46 l 
‘ . 9*? 49 4 
) . 99 V^4 


C .8 

r. 45455 0 
,.46635 0 
f .47 2 24 C 
0.47812 0 
0.43398 0 
2. 4 8 9*2 > 

c. 40564 0 
r . 50 1 44 C 
7. 51294 C 
2.52432 0 
C. 53 55 7 0 
j . 64 6 7 ? C 
C. 55 766 0 
0.55341 0 
0.615'' 7 0 
0.66413 C 
*.70 36 0 9 
2 . 7 ? 7 o 7 0 
2.74697 C 
2.79213 0 

? . 8 3 2 3 6 0 

0.36249 0 

:.« { *541 3 

° 30 9 0 

C . 9 1 7 11 C: 
0.92933 C 

* , o 4 r 7 7 

r.97534 : 

} . 5 ° 3C H 0 
'*.93 76? 0 
0.99060 3 
0.99261 0 
■*.99 50 9 : 

: . o 9 7 3 8 : 
0.99A67 0 

:. 99965 c 

* . 9 > J94 '* 

0.99904 r 


C.7 
43 73 l 

5 : c 1 0 

,506 21 
, 6 1 ? 7 9 
,6 1334 
,52435 
,5303 3 
,53627 
,5490 ) 

, 55955 
, 57 C9 l 
. 59206 
. 59299 
.6^373 
.65341 
.4947? 
.73391 
. 7*561 
.79256 
.81550 
. 85134 
.97879 
. 39 9 1 C 
.91436 
.92 724 
.9371? 
.96564 
. 97841 
.93511 
.83917 
.99177 

. 99354 

.99570 
.99771 
. 99384 
. 99972 
. 9099 5 
.99995 


0.6 

C .5 *>63 2 
0.539C2 0. 
3.54530 ‘ 
j. 55 154 0. 


_ 58 772 
0.56385 0. 
;.5o992 
C . 5759 3 0 
0.5 377 5 0 
C.4oqi3 C 
0 • 4 1 C6 4 0 
C. 42167 J 
0.63242 0 
2.64233 0 
:. 49081 0 
0 . 73149 0 
w. 76597 0 
0 . 79435 0 

C. 31 91? 0 
0.83959 0 
0.K7173 ; 

0.39535 ■: 

C . 9 1 3 0 ft 0 
0 . 92 o 7 3 ■* 

. 93 74? O' 
0 . 2] 4 6 5 7 0 
C. 97C5 2 0 
0.58149 O 
0.93724 
0.9907? 

J. 99 2 95 
0.99446 
0 .996 3? 

0 . *>9 3 04 
0. 9990 1 
j . 99974 
C .99996 
9.9999a 


0.5 

- 57143 
53433 
59075 
597C5 
6 0 3 2 4 

60 9 42 

61 546 
67144 
6 3 314 
7)4 45 1 
45654 

« 5 5 6?? 

. 5766s 
.43654 
.73158 
. 76396 
.79934 
. 8? 5 4 7 
. 34467 
.36 4 ? 6 
. 847C4 
.912 U 

. 92 7 1 7 
. o 3 3 6 9 
, 94 7 79 
.95436 
. 97541 

. •76<»5 7 

.98 n 7 

.99777 
. 9 74 1 3 

, 94537 

. 49405 
. 94 8 3 6 
. Q <J o 1 7 
. 94978 
. 49496 
. 999 9,4 


0.4 

.62500 
.63790 
.6R4 2C 
.65079 
.65643 
.66247 
.6h335 
.67412 
,s3535 
.69616 

. 70 6-57 
.71456 
. 776 1 7 
. 7 3 5? 5 
.7761 * 
.8" 70 L 

,89467 
.847*0 
. P 7 5 1 4 
.889 9 6 
•°1273 
.9?619 
.541 4 3 
, 091', 7 6 
. <347:2 
* 9 6 7 7 0 
1 . 9 3 0 '* 1 
1.93 75 7 
. 9 9 l 6 C 

‘ . 993 M ? 

. S453C 
), 994 3 1 
>.99755 
). 99369 
'.99974 
; . 9498 3 

■ .00947 

99999 


0.3 

0 . o 8 96 6 
0.70195 
0.70790 
0.71371 
0.71939 
C. 72495 
C . 73037 
' .7 3566 
2 . 745 Q 3 
t . 75 5o 2 
0.76439 
C . 7737? 
C . 7 3 2 1 3 
0 . 79C 13 
0.32475 
C . a 5 l p 4 
0.37343 
0 . 5 9 C 7 4 
0.90479 

C . c 1 6 4 2 
C . 9 3 1 9 2 
; .9444 3 

6 Q 5 ^ ^ 4 
? „ 5C?< 0 
0 .96842 
r . o 7 q 

7 .5 3 >22 

0.49069 

C . o q 1 , * 3 
C . 9 9 5 3 7 
0/ . 9 9 r> 4 6 
0.99773 
<’.99816 
0.999 0? 
0 . 95 95C 
0 . 99 73 7 
<*.99973 
j. 99^48 


0.2 

C • 7 69 ? 3 
0.77934 
0. 7H4Q7 
C. 7399C 
0.79465 
0.79934 
C. 80393 
0.80831 
0.81648 
0. 8?456 
2 . » 3157 
C .« 38 7 5 
0.84553 
C.85L74 
0.37 794 

0.89791 
0.51341 
2 ,9?S4'i 
0 . 5 354 6 
0.94345 
0 .9556 3 
r . 96‘t 1 Oil 
2, . 9 7 0 4 1 
0.3751'' 
0.97862 
0.581 77 
C. 99007 
# q O 3 3 ^ 
0.94576 
0.99691 
0.99765 
■3,95816 
0.94877 
9.99935 
C .999b 7 
0.99971 
0.99991 
0.54991 


0. I 

0 . 86957 
0.87671 
0. 8«C09 
C . 98332 
0 . 88643 
C . 8894 l 
0.85228 
C. 99504 
0.90021 
C.9C504 
C. 90452 
0.91 169 
0.91758 
0.92121 
0.93621 
0.94722 
0.95549 

0.96204 
0.96710 
0.97127 
b. 07747 
0.98188 
0. 98 5 1C 
C. 98752 
* .909 I® 
0.99087 
C . 99 5 05 
j .9 46 90 
0.90788 
0.99846 
0.99881 
C . 99938 
C . 999 19 
0.99967 
0.99904 
0 . 99996 
3 . 99996 
0.99996 
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